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CHAPTER I 


INTRODUCTION 

The Research Grant NGR-06-002-147 is concerned with the system 
analysis of collision-dominated and collisionless plasma centrifuges 
and the theory of sputtering and deposition of sputtering products on 
surfaces. The period of the grant extended from 1 November 1976 to 
31 October 1977. This is the final report on the research carried 
through in this period. 

In Chapter II, an analytical theory is developed describing the 
deposition of sputtered atoms on system surfaces which cannot be seen 
along straight paths from the emitting surface. The boundary-value 
problem describing the diffusion of the sputtered atoms through the 
surrounding rarefied electron-ion plasma to the "hidden" system surfaces 
is formulated and treated analytically. It is shown that outer boundary- 
value problems of this type lead to a Fredholm integral equation. The 
latter is solved by the method of successive approximations. A quantum 
theory of sputtering of metal surfaces by low energy ions (<100eV) has 
been developed and submitted for publication. This work will be commun- 
icated at a later date. 

In Chapters III and IV, centrifuge models employing ring electrodes 
of different radii located in the end plates of a cylindrical discharge 
chamber, are analyzed which avoid the boundary layers at the inner electrode 
cylinder and (probably) the secondary flows and instabilities occurring in the 
magnetohydrodjmamic flow between concentric cylinders (Chandrasekhar 
1961). Complete two-dimensional solutions are derived which show that 
the Hall effect enhances plasma rotation (Chapter III) and that the 
Induced magnetic field does not interfere with the rotation due to the 
external magnetic field (Chapter IV) . These schemes exhibit velocity 
end losses due to the boundary layers in the cathode and anode planes. 
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For this reason^ an. Improved centrifuge system is conceived, C7hlch 
essentially avoids the velocity end losses caused by boundary layers 
at the electrode plates (Chapter V) . In view of the circumferential 
electrode arrangement, a multidischarge counter-rotating centrifuge can 
be set up in a long insulating cylinder to rotate large volume of isotope 
mixtures. In Chapter VI, the theory is applied to the separation of 
^238 ^235* difficult problem of compressible plasma centrifuge 

analysis is formulated in Chapter VII and solved by means of Lyapunov- 
Schmidt series expansions. In Chapter VII, a simple theory for a colli- 
sionless plasma centrifuge is formulated based on the coupled Vlasov- 
Maxwell equations for the electron and ion components. 
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CHAPTER II 

SOLUTION OF EXTERNAL BOUNDARY-VALUE PROBLEM 
FOR DEPOSITION OF SPUTTERING PRODUCTS 

In an ideal vacuum* sputtered atoms travel undeflected along straight 
paths determined by their initial velocities at the point of emission. 

Within this free particle flow* a system surface is reached by the sputtered 
atoms only if it can be seen along a straight line from the emitting sur- 
face. In reality* ion propulsion systems are surrounded by a very rarefied 
plasma consisting of escaped beam ions* recombined ions, and electrons. 

For this reason* always some of the sputtered atoms will be deflected out 
of their initial paths by Interacting through long-range forces (polariza- 
tion forces) with the plasma particles so that they can reach system sur- 
faces which are not seen along a straight line from the emitter. 

An idealized propulsion system exhibits an emitting plane 

^ ® ,1 ^ ^ (accelerating grid), the rocket surfaces r ■= a, -c £ z £ 0 

and z = -c, 0 £ r £ a, and the plane z = -d* a^rj<bof the solar energy 
collectors. All these system surfaces can be reached by the atoms sputtered 
from the emitter by diffusion through the rarefied plasma. The diffusion 
coefficient D is determined by the Vlasov equation for the sputtered atoms 
Interacting through weak long-range forces with the plasma particles. In 
view of the mathematical difficulties associated with the solution of 
boundary-value problems for this geometry, a somewhat simpler system is 
studied here consisting of an emitting plane (z = 0, 0 r a) , the upper 
rocket surface (r = a* -c £ z ^ 0) and the plane (z “ -c* a ^ r <^ «>) of 
the solar energy collectors (Fig, 1). The latter is assumed to have 
infinite radial extension. 

Ordinary boundary-value problems are defined for a space bounded 


on the "outside" by boundaries, whereas external boundary-value problems 
are defined for the space surrounding "inner" boundaries. In terms of 
cylindrical coordinates (r,0,z), the space in Fig. 1 consists of the 
adjacent regions, 

I; 0 z 0 r 
II: -c ^ z £ 0, a ^ 

In this case, the "inner" boundaries are formed by the cylindrical wall 
(r = a, -c _< z £ 0 ) and the circular end-surface (0 _< r £ a, z = 0 ) of 
a cylinder of radius r = a extending from the plane z = -c to the height 
z “ 0. On the other hand, the plane z = -c, a^r^“ represents an 
external boundary of the space II. 

We consider herein the external boundary-value problem for the steady- 
state diffusion (Laplace) equation and the space I + II shown in Fig. 1 
when the end-surface (0 ^ r ^ a, z = 0 ) of the cylinder emits particles 
at a given rate I(r). At the inner, cylindrical boundary (r = a, 

~c £ z £ 0 ) and the bottom plane (z * -c, a £ r ^ ®) the particles are 
assumed to be deposited by adsorption or absorption. 

Various other transport processes for particles or heat in technical, 
physical, and biological systems lead to external boundary-value problems 
of this type. We mention as examples 1) the emission of particles from 
a cylindrical chimney into a gaseous atmosphere, and il) the Injection 
of a liquid from a cylindrical probe into a biological medium. 

Analytical or numerical solutions of external boundary-value problems 
have apparently not been given in the literature. We will demonstrate 
that the considered external boundary-value problem can be solved analyt- 
ically by means of a Weber transform. In this analytical solution a 
matching function i()(r), a 5 r js « (at the interface of the regions I and 


i'..- V. 
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II), occurs which is determined by an inhomogeneous Fredholm integral 
equation of the first kind. This integral equation is discussed and 
transformed into an inhomogeneous Fredholm equation of the second kind, 
which is solved by the method of successive approximations.—^ 
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In the space z ^ -c, let the density of the diffusing particles be 

-3 

designated by n(r,z) [cm ] and the flux of emitted atoms at the emitter 

-3 —1 

surface by I(r) [cm *cm sec ]. In steady state, the spatial distri- 
bution n = n(r,z) of particles is determined by the external boundary- 
value problem for the Laplace equation (Fig. 1): 


[9n(r,z)/3z]^^g = -I(r)D 


0 < r < a 


n(r,z) = 0, -c < z < 0 
r=a — 


n(r,z) = 0, a < r < “ 
z=-c * — — 


2 2 

n(r,z) 0, (r + z ) ’>■ 


as the proper and improper boundary conditions, respectively. D 
designates the diffusion coefficient of the particles. 

The boundary conditions (3)- (4) imply that particles arriving at 
the indicated surfaces are deposited there, l.e., do not return into 


|!1 


• V 

N _ = -2iraD / [3n(r = a,z)/3rj( 


BOUNDARY-VALUE PROBLEM 


i'i- ! 


a^n n 

. 2 r 3r 2 ^ 

dr 3z 


the diffusion space. The fluxes = -D V^n of particles arriving at the 
system surfaces r = a, -c £ z £ 0 and z = -c, a£r_<» are given by 
(Fig. 1): 


1 I 

. f . I 


= a,z) = “D 3n(r = a»z)/3r, “C £ 2 ^ 0 

^ (2 " -c,r) = “D 3n(z - -c,r)/3z, a < r < ® 

z ~ 


Aci'ortllii}:iIy, 
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• « 

Af __ = -2irD / [3n(r,z = -c)/9z]rdr , (9) 

a 

are the numbers of particles deposited per unit time on the system 
surfaces r = a, -c j< z ^ 0 and z = -c, a _< r _< respectively. 

The above boundary- value problem cannot be solved directly, i.e., 
requires a decomposition of the space z ^ -c into appropriate "subregions 
I and IT for which the associated boundary-value problems are solvable. 

In this approach, a common boundary value [i/»(r)] at the decomposition 
plane is determined by an integral equation. 

Let dimensionless independent and dependent variables be introduced 
in accordance with: 

p = r/a, 0 _< p _< ? = z/c, -l ^ 5 £ m , (10) 

and 

N(p,?) = n(r,z)/n , S(p) = I(r)/I , (11) 

o o 

with 


n = cl /D, I = l(r=0), y - c/a 
o o o ' / » 

The boundary-value problem defined in Eqs. (l)-(5) reads in 
dimensionless form: 


( 12 ) 


9"N . 1 3N . -2 9^N „ 

77‘^p9^""Y 7T=° 

9p ^ 3? 


(13) 


where 


[3N(p.?)/3?]^^q = -S(p), 0 < p < 1 , (14) 

N(p,5)p^j^= 0, -1 1 C 1 0 , (15) 

N(p,C)^^_j^ = 0, 1 < p < ® , (16) 

and 

N(p,C) 0, (p^ + 5^) » . (17) 


In Fig. 1, the space is decomposed into the regions 




I 1 I 


1 ^ v> 

\1 ■: 


1(0 ^P£“» 0^5_<<») and 11(1 p ^ -1 £ ? £ 0) . At the interface 

? = 0, 1 £ p < «, the partial 3N(p,C = 0)/3? = ¥(p)H(p-l) is introduced 
as the conmion (unknown) boundary value 4'(p) of the adjacent regions I and 
II, 1 £ p £ “, Thus, the boundary-value problem in Eqs. (13)-(17) can 
be decomposed into boundary-value problems for the regions I and II. 

I. In region I, N = Nj.(p,?) is described by the ordinary boundary- 
value problem; 



-2 

+Y _=0, 0<p<“, 

3? 

9 

(IB) 

[3Nj(p,0/9cI 

5=0 " 0 £ p £ 1 - 0 

f 



= ’f’(p) » 1 + 0 £ p £ 

9 

(19) 

N^(p>?) 0, 

(p + 5 ) + » 

9 

(20) 

S(p>l+0) = 0, 

^(pcl-O) H 0 


(21) 


where 


for physical reasons. Since region I is the upper half of the infinite 
space (0 < 5 £w), the general solution of Eqs. (18) and (20) is given 
by the Fourier integral, 

03 

N^Cp.C) = / A(k)e~^^^ J (kp)dk , (2: 

0 ° 

which satisfies the improper boundary condition for p -»- » and ? 

The Fourier amplitude A(k) is determined by the boundary condition (19) , 

CO 

-Y / A(k) J (kp) kdk = -S(p)H(l-p) + P(p)H(p-l) . (2: 

0 ° 

Application of the inverse Hankel transform to Eq. (23) gives 
1 -1 ® 

A(k) = Y / S(a)J (ka) ada - y / 'l'(a)J (ka) ado . (2i 

0 ° 1 ° 


r 
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Substitution of Eq. (24) into Eq. (22) results in the solution for 
region 1: 

00 1 eo 

N^(P,5) = / dk J (kp)[/ S(a)J (ko) ada - / 'P(a)J (ka) ada] , 

^ 0 ° 0 1 

II. In region II, N = Nj.j(p,?) is described by the external 
boundary-value problem: 


(25) 


+ i !!!ii + ° "ii , 0 

(3N3-j(p.O/3C1^^0 = 'J’(P), 




1 ^ p ^ “ , 

-1 5 C 1 0 , 

1 ^ p ^ , 

p -V «», ”i ^ C ^ 0 

According to Eq. (28), region II has an inner, cylindrical boundary 


^Il(P,^)p=l 

Nij(p,i^)^^_l - 0, 

Njj(p,0 0, 


(26) 

(27) 

(28) 

(29) 

(30) 


at p = 1 where N^j(p,c) vanishes. For this reason, a Fourier integral 
representation of Njj(p,c) is needed for 1 ^ p ^ « which vanishes at 
p = 1. According to Weber’s integral theorem, an arbitrary function 

2 3) 

f(p)» ft P < ®, with f (p = ft, “) = 0 satisfies the integral equation:—’— 


® kdk Wj'(p) “ 

f (P) = I -o ^ / f(ot)W^(o) ada, ft 1 P 5 

0 J^(kft)+Y^(kft) ft 


(31) 


where 


wj;(p) H J^(kp)Y^(kft) - J^(kft)Y^(kp) , (32) 

and J (kp) and Y (kp) are Bessel functions of order v of the first and 

V V 

second kind, respectively. In view of Eqs. (31)-(32), a Fourier integral 
solution of Eqs. (26)-(3Q) is sought in the form. 



r!”T Tin 


FH Ml 1 II 




^ * 
i i 


Cp) sitih[Yk(i;+l)ldk 


Wj.(p) 5 Jjj(kp)Y^(k) - J^Ck)Y^(kp) 


which obviously satisfies the boundary conditions (28)- (30). The 
Fourier amplitude B(k) is determined by the boundary condition (27) , 


Y / B(k)W. (p' coshYk kdk = ^(p), 1 P 5. " 


which gives 


B(k) = {ycoshYk* [J^(k)+Y^ (k) n ^ /*if(c!)W,.(«) ada 


o o 


by Eq, (31). Substitution of Eq. (36) into Eq. (33) results in the 


solution for region II; 

TJ fn r'i f n- sinh[Yk(?+i)] . 7 

II Q cosliYk j j 


/ ¥(a)W. (ct) ada, 


l£p£co, 

The solutions N (p,?), Eq. (25), and Nt^(p,^) Eq. (37), contain 
I 

the yet unknown boundary-value ¥(a), 1 ^ a £ <». ¥(a) is determined 

by the continuity condition at the interface of regions 1 and II, 
Nj(p,5)^„g = 1 1 P 1 ® 

which gives 


/ dk J (kp)[/s(a)J (ka) ada - / ¥(a)J^(ka) ada] 

0 o 0 ® 1 ° 

® W (p) w 

= /dk tghYk —5 5 /’P(a)W- (a) ada, 1 £ P i “• (^5) 

0 J^(k)+Y^(k) 1 

Eq- (39) indicates that '?(a) is determined by an inhomogeneous integral 

1 ) 

equatxon,— 




. ...t.:-: I . 1... 




‘ ^ • 


Because of the boundary conditions (19) and (27), the remaining 
continuity condition at the interface of regions I and II, 

[3Nj(p,5)/3?1^^q = lSNj.j(p,?)/a?]^^Q, 3-lPl" . C 

has already been satl led. Indeed, substitution of Eqs. (25) and (37) 
into Eq, (40) yields 


“/ kdk J (kp)[/s(a)J (ka) ada - /v(a)J (kc:) ada] 

0 0 ^ 1 ^ 


= f kdk 


,2 2 , 


0 r(k)+Yj(k) 1 


/ 'y(a)W. (a) ada, 1 _< p 


By means of the Hankel and Weber integral representations for the Dirac 
function 6 (a-p), 


/j^(fcp) J^(ka) kdk = 6 (a-p)/a 


“ W (p)W. (a) 

/ —5 —5 kdk = 6 (a-p )/a 

0 J^(k)+Y^(k) 
o o 


Eq, (41) is reduced to 


- /s (a) 6 (a-p) da + /'i'(a) 6 (o-p) da 
0 1 


= /f(a) 6 (a-p) da, 1 £ P £ ® 

1 

Eq, (44) gives the expected identity, ’i'(p) = V(p), since the first 
integral vanishes by Eq. (21). 


JF-- 


INTEGRAL EQUATION 




By introducing the kernel fCCa^p) and the source Q(p), the integral 
equation in Eq. (39) can be revrcitten in the convenient form: 

CO 

/'P(o)K(a,p) ada = Q(p), 1 ^ P 5 ” » CA5) 

1 

where 

“ « W. (a)W (p) 

fC(a,p) = /j (ko)J (kp)dk. + /tghvk-s 5 dk , (A 6 ) 

0 ° ° 0 J^(k)+Y^(k) 

«o 1 

Q(p) = /dk J (kp)/s(a)J (kcc) ada . (47) 

0 ° 0 . ° 

For simple particle emission distributions S(a), e.g., in the case of 
a homogeneous and a parabolic emission distributions, respectively, the 
source integral Q(p) is readily evaluated, 

Q(P) = f P [E(-^) - (1 - -”)K(^)1 , 1 1 P 1 “ 

for S(a) = 1, 0 a _< 1 , (48) 

and 

Q(p) = ^ [(4-2p^)E(|)-( 1- -^)(3-2p^)K(|)] 

for S(a) = I -a^, 0 ^ i ^ 

1 1 

K(“") and E(— ) are the complete elliptic integrals of the first and 
P P 

second kind, respectively [K(l) = E(l) = 1].—^ Eq. (45) reduces 
the external boundary-value problem to the resolution of an inhomogeneous 
Fredholm integral equation of the first kind for the unknown boundary- 
value V(p), 1 ^ P £ 

Comparison of the integrals in Eq. (46) with Eqs. (42) -(43) 
indicates that the kernel K(a,p) is the sum of two integral functionals 


zr 


■" T- 
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which are singular at a *= p, similar Co the Dirac function. It follows 
that the solution of the integral equation (45) is ¥(p) = 'PqCp) = j Q(P) 
in the lowest approximation. By means of Eqs, (42)-(43), we transform 
Eq. (45) into the Fredholm integral equation of the second kind:— ^ 

CO 

’i'(p) = Q(p) jH'(a) K(ct,p) da, 1 < p £« , (50) 

^ 1 


where 

, « “ W, (a)W. (p) 

K(a,p) = a/j (ka) J (kp) (k-l)dk + a / — = x (k-cghyk) dk , (51) 

Q° ° 0 J^(k)+Y'^(k) 

c o 

Q(p) = j Q(p) • 

It is seen that the parameter of the integral equation (50) is X = 1/2. 
The kernel K(a,p) in Eq. (51) consists of two integrals, each of which 
is the difference of two integral functionals which go to “ for a + p. 

Eq. (SO) is solved by the method of successive approximations which 

.1) 

gxves:— 

'f(p) ** lim l'jj(p) » 1 1 P 1 “ » 

n-K» 

where 

4*q(p) - QCp) 


'I'j^Oj) = 'i’^Cp) + c|)^ / Q(a^)K(aj^,p)dUj^ 

*^Tlie singular behavior of K(a,P) is demonstrated by evaluating the 
integrals in Eq. (46), e.g., 

“ -1 

/j (ka)J (kp)dk = (2/tt)p K(a/p), 0 < a < p < ® 

0 ° ° 

= (2/Tr)a'"^K(p/a), 0 < p < a < «» 

where 

K(l) « » 


” ^ ^ J - J 1 . j . 1 - — > 1 — - j — j 
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1 1 


00 CO 


1'^(p) = 'l'n_i(p) + (•!)"/... /Q(ct^)^(a^,a^_j^)...K(ct2,otj^)K(aj^,p)da^...daj^ 

(54) 

Combining of Eqs. (53) and (54) yields the solution of the nth 
approximation in the form; 


n 


(p) “ 2 U, (p)y tt“ 0>lj2)»«« ™ J ^^P ^ 

^ i=0 ^ 


(55) 


where 


Uq(p) = Q(p) 


Ul(p) = (|■)^/p^(aJ^)K(aJj^,p)da^ 


U2(P) “ 



00 00 ^ ^ ^ 

/ / Q(a,)K(a^,OL)K(a 

11 " ^ i i 


,p)da^dQj^ 


f 


Uj^(p) = (■|)'^ /•••/ Q(a^)f((a^,ajj_j^) • . .K(a2,a^)K(a^,p)dtt^. . .doj^ 


(56) 


It is readily shown that the Neumann series in Eq. (55) converges 
to 'l'(p) in the limit n •+ ™.— ^ It should be noted that *4' (1+0) is not 
necessarily S(l-O) in any approximation n, i.e. 3N(p,0)/3c may be dis- 
continuous at p = 1. With ^Cp) given by Eqs. (55)-(56), the particle !/• 

/ 

density fields N^(p,?) and Nj^(p,c) are given by Eqs. (25) and (37), 
respectively.. Thus, we complete the analytical solution of the external 
boundary-value problem under consideration. 
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CHAPTER III 

PLASMA ROTATION BY LOKENTZ FORCES 
IN A DIVERGING PLASMA CENTRIFUGE WITH HALL EFFECT 

In this chapter, a system analysis for a collision-dominated 
plasma centrifuge is presented in which the plasma rotates under the 
influence of the Lorentz forces due to the interaction of a spatially 
diverging current density field with an axial external magnetic field. 

The associated boundary-value problem for the coupled partial differ- 
ential equations, which describe the electric potential and the plasma 
velocity fields, is solved in closed form. The electric field, current 
density, and velocity distributions are discussed in terms of the 

Hartmann number H and the Hall coefficient ut. As a result of the 

4 4 , 

Lorentz forces, the plasma rotates with speeds as high as 10 m/sec 

around its axis of symmetry at sufficiently large values of H and ojt. 

It is remarkable that the Hall effect supports the plasma rotation. 

3.1. Model of Diverging Plasma Centrifuge 

An electrical discharge in a cylindrical container rotates if the 
Lorentz force has a nonvanishing component in the azimuthal direction. 

For example, arc experiments in an axial external magnetic field B 

2S 

indicate that the discharge plasma rotates (Schwann 1970; Vedenov et al. 

1961) since the current field lines J have a nonvanishing radial 

component J so that (JxB)q=-JB ?^0. Ina stable arc discharge, 
r 0 r 2 


. •. 





1 i. X / . i ' - i. - -X i 


the J^-component Is caused by the concentiration. of field lines (J) 
the electrodes and a dilatation (repulsion of currents in the same 
direction) of the field lines J in the Interelectrode space. 

A theoretical model for the production of a high-density plasma 

centrifuge, which has a radial current density which is In magnitude 

comparable with the axial current density J , is shown in Figure 3.1. 

The radial spreading of the current field lines J is forced by means of 

electrodes of considerably different radii R, and R (R. » R ) in the 

end plates z ** ±c of an electrically isolating centrifuge chamber of 

radius R^. The field line? of the current density J and of the external 

axial magnetic field cross under a nonvanishing angle (except at the 

chamber axis) so that the resultant Lorentz force J x B rotates 
^ o 

charged particles around its axis of symmetry. In the steady state, 
the magnetic body forces in the azimuthal direction are balanced by the 
viscous forces (boundary layers at the chamber walls) . As opposed to 
the centrifuge with radial electric current flow between an inner and 
outer cylinder electrode, the centrifuge scheme in Figure 3.1 avoids 
the boundary layer and losses at the inner cylinder surface. 

In the following, the steady-state rotation of the spatially 
diverging plasma contained by an insulating cylinder in the external 
axial magnetic field (Fig. 3.1) is treated theoretically. The 
analysis is based on the magnetogasdynamic approximation, in which two 
characteristic dimensionless parameters occur, the Hartman number H 


and the Hall coefficient oix. The magnetic fields associated with the 

discharge currents « {J , J , J } are neglected for small magnetic 

H b Z 

Reynolds numbers [a, = o/ (l-ha^T^) ] , 




i. I I Mil ! 
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<< 1 , 

Rg « 0[Bg/B^I » lu^I/RirR^j/B^ « 1 , (3.1) 

R « 0[B /B 1 = n WTO.V R « 1 , 
z *zo o-**oo 

« 

vhere is the characteristic velocity of rotation and 1 the discharge 

current. These inequalities are satisfied in many cases, e.g., if 

R 0 « Ifor i) WT » 1 and R < 1 or ii) R « 1 and 0 ^ wt < «»^ where 

R s u ov max (R ;c) , 
o o 

3.2, Boundary-value Problem for Velocity and Electric Potential 


For a purely azimuthal flow field V » {0, Vg(r,z), 0}, the plasma 

behaves incompressible, V * V * 0. From the continuity equation (2.2) 

for the steady state, V *(p V) = V • Vp =» 0, it follows then that the 

P P 

density gradient Vp^ is everywhere perpendicular to the flow field V. 

These ideal conditions are realized if secondary flows are absent or at 

least negligible (Schlichting 1960). In accordance with the steady-state 

magnetogasdynamlc equations (Eqs. (2.1), (2.7) and (2.8)], the rotating 

plasma in a homogeneous magnetic field (Fig. 3.1) is described by 

the following boundary-value problem for the azimuthal velocity V.(r,z) 

□ 

and electric potential fields [induced magnetic fields neglected, 

Eq. (3.1)]: 


a 1 a 

8 r ^r ^ ^ 

1 _B_ / ^ 

r 8 r 30-' 


(. M + V B 1 

B Mt ^ Vo' > 


V 


(3.2) 


3 


(3.3) 




. i . .. . I . - .1 


where 


Vg(r,2:)^^j^ = 0 , -c £2 < + c, 

o 

= 0 * 0 1 :i ^o * 


(3.5) 


[3(f(r,z)/3r]^_j^ =0 , -c <_ z <_ + c. 


(3,6) 


ff(3(J(r,z)/3z] _ . * I 6(r - R )/27rr , 0 £ R • (3.7) 


The boundary conditions (3.4), (3.5) and (3.6) consider that the plasma 
does not slip at the walls r ~ and z = ±c, and that no current flows 
into the cylinder wall r ” R^, respectively. The boundary conditions 
in Equation (3.7) imply that the cathode (R_) and anode (R_j_) are ring 
electrodes of vanishing radial width, Ar -»■ 0[6(r - R^)/2in: = radial 
Dirac function] . The net current flowing through the centrifuge is by 
Equation (3.7). 

-2to f ** 9<KryZ=±c) . ^ r o dr = I < 0 , (3.8) 

JQ dZ ■'Q ± 


since the positive current (X < 0) flows from the anode to the cathode 
(Fig. 3.1). The pressure distribution P *» P(r,z) is determined by the 
r^* and z-components of the equation of motion. 


“p “ = - ~ + aiT o, B (- 1^ + V B ) , 
p r 3r -*• o 3r 0 o 


(3.9a) 


Ve- Vr-^“* 


(3-9b) 


According to Equation (3,9b), it is 3P/3z = 0 if the induced fields B^ 

and B_ are neglected [Eq. (3.1)]. This means that momentum cannot be 
0 




exactly balanced In the axial direction If Induced magnetic fields are 

neglected (in absence of secondary flows). Equation (3.9b) is in 

accord with the boundary-layer approximation according to which the 

normal pressure gradient is 3P/3z - 0 at the electrode plates s = ±c, 

0 < r < R . 

— — o 

In the absence of the Hall effect, iot « 1, it is V x B = u 

o 

(J , 0, J ). Hence, B =0 and B =0 because of the homogeneity of 
r z r z o J 

the boundary conditions for B and B , whereas B 0 since J # 

1C Z O TC ^Z 

0[Bo(r,z=±c) = tp I/2Trr)W(r-R.)l . Consideration of the Induced field 
□ Ox 

B = {0, Bg, 0} leaves Equations (3,2) and (3.3) unchanged. This means 
that the boundary value problem in Equations (3.2)-(3.7) and the 
solutions Vg(r,z) and ^(r,z) derived from it remain valid even in 
presence of a significant induced field B = (0, Bg, 0}, Rg i 1, as long 
as the Hall effect is negligible, ut « 1 (Chapter IV). 

3.3. Fourier-Bessel and Dini Series Solutions 


The characteristic dimensionless parameters of the magnetogas- 
dynamic centrifuge problem under consideration are obtained by 
intorducing the dimensionless independent and dependent variables. 


p « r/R^ , 0 p £ 1, (3.10) 

C « z/c , -1 1 ? 1 1» (3.11) 

and 

V(P.C) » Vg(r,z)/V^, *(p,C) = (3.12) 

where 


= ^^/R B , « 5 Ic/2iictR^ . 


(3.13) 


In terms of the dimensionless space variables and fields, the 
boundary-value problem defined in Equations (3. 2) -(3. 7) assumes for 
V(p,0 ^nd $(p,z) the form: 




a« 


,-2 3^ _ _ _g2 

2 1 1 3p 


(3,14) 


35“ 


kJL 4. m~ 2 3^ - i A. 

p Dp 3p^ g^2 p 3p ^ 


(3.15) 


where 


and 


with p^ 


V(p,5)p_j^ « 0, -1 < ? < + 1 , 

V(P,?)^^ “ 0, 0<p<l, 

t3®(p,C)/3plp^i = 0 , -1 1 Cl 1 » 

-[3®(p>C)/ac]^=^ = 5(p-P+)/P , 0 IP 11 . 

R./R . The dimensionless constants M, N, and H 
± • o 


(3 

(3 


,16) 


.17) 


(3.18) 
(3.19^ 
are defined 


by 

M 5 (1 + (c/Rp), N e c/R^ , 

H s (a /u)^ 3R ^ H /(1 + uj2t:2)^ 
11 ° ° 



In view of the similarity of the left sides of Equations (3.14) 

and (3.15) with Bessel's differential equation, Z" + p ^E' •{• 

m m 

9 “2 2 

(fc2 - p m^)Zj^ = 0, for cylinder functions ^^^^(k^p), partial solutions 
of the coupled inhomogeneous equations are sought in the form. 


V^(p.C) - Ji(V^^n^^^ ’ 


*„(P,C) - Jntk„P)8„(C) 


(3.21 




(3.22 


u 

where Ji(k P) » -J. (k p) and J*(k p) + (k p)“^J. (k p) » J,,(k p). 

On In in n in Un 

Substitution of Equations (3.21) and (3.22) into Equations (3.14) and 

(3 .13) yields 

f” - + h^)n\ = k^^aVg^ . (3.23) 

^ - k^ . (3.24) 

where the eigen-values k^ > 0 are detennined by the boundary conditions 
(3*16) and (3.18) as the real roots of the transcendental equation 

Jl(k^) = 0 , Si = 1,2,3,.... (3.25) 

Thus, the general solution of the coupled equations (3.14) and (3.15) 
obtains by linear superposition as 

OD 

V(p,c) = S Ji(\p)f„CC) . (3.26) 

n=l 

CO 

«(P,5) = -2? + E J (k p)g (CJ • (3.27) 

n-1 ° ^ 

In view of Equation (3.25), Equation (3.26) is a Fourier-Bessel series, 
whereas Equation (3,27) is a Fourier-Dini series in which a zero-order 
tema, -2^ has to be included, in accordance with the Fourier-Dini 
expansion [Eq. (3.32)] of the boundary value in Equation (3.19). By 
decoupling Equations (3.23) and (3.24) one finds for and g^(c) 

the differential equations of 4th order, 

f’"’ - [k^(M^ + N^) + N^H^Jf" + k^M^N^ f - 0 , (3.28) 

n n n n n 

g"" - tk^(M^ + N^) + NV]g” + g * 0 , 

n n i n “ 


(3,29) 
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t i. 1- L 


with 






C3.30) 

C3.31) 


as boundary conditions by Equations (3.17) and (3.19). In deriving 
Equation (3.31), the Dirac function in Equation (3,19) has been 
-expanded as the Fourier-Dini series. 


«<P P+)/p -2 + 2 S [JqOc p^)/J^(y]J (k^p). 

n-1 


1(3.32) 


In addition to Equations (3,28)-(3.31) , sund S^^C?) have to satisfy 

also the coupled Equations (3.23) and (3,24). With 


+ — + 

0- *= (»3 , a- = u , u„ = -u , w. 

In n * 2n n * 3n n* 4n 


n 


(3.33) 


where- 


“n = \j + N^) + N^H^] ± {[k^(M^ + N^) + 

- 4kViI^r'}]^ , (3.34) 

± + 

W C Ci)~? 

the general solutions for f^(?)“ e ^ and g^(C)« e of Equations 
(3.28) and (3.29), can be written as 


"f* 

. sinh w ? . cosh oi ? 

f (5) « I-+B'*' ^ 

n n .+ n + 


sihh to 


cosh to 


(3.35) 


+ A 


sinh to C 

rr 


sinh to 


cosh e Z, 

+ B*“ ^ 

n 


n 


cosh 


^ 9 


±__,L_L.. 


ilA 


.1.. 


i X. 
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•f* *f* 

. s±cih (d C , cosh w C 

^ f- 

slnh (d cosh u 

n n 


+ C 


slnh £d~C 

— n 


n 


+ D 


sinh (d 


n 


cosh td C 


cosh (d 


(3.36; 


r ! 


Only four of the eight integration constants D“ for any n > 1 

n n — 

are Independent; by Equations (3.23) and (3.24), 


and 


[Ctd“)2 - - k , 

Ti -rt n 1-1 ^ 


n n n 


[(w*;2 - k^M^]D“ = fc , 

n n n n u * 


(3.37) 


[(/)2 

n 


u 


(k^ + H^)N^]A“ = k , 

n jL n n j.n-* 


Ck“ + H^)N^]B" = k K‘-HrD" , 


2 

n 


2 2 -f- 
rE D" 
n ^ a 


(3,38) 


where the coefficient determinants of the pairs of corresponding 
equations in Equations (3.37) and (3.38) vanish owing to Equations (3.33) 
and (? 34). 

«*• 

Upon application of the four relations in Equation (3.38), which 
•are equivalent to Equation (3.37) by Equations (3.33) and (3.34), and 
the boundary conditions (3.30), which give 


-a" * +A‘‘* s a , 





(3.39) 


Equations (3.35) and (3.36) become 




C slnh Id C sinh id"c\ 

J 

slnh (0 sinh u / 


+ B 


i 


cosh (d 5 cosh 03 z 


n 


nv 1. + 
cosh u 

n 


cosh 03 


r) 

n 




n 


k 

n 1 


B 


n 


k 

^ 1 


{< 

{< 


sinh os"*"? 


sinh 0 ) 
n 


“f- 

cosh 01 Z 
^ 

cosh 03 ”*" 


sinh 0 ) Z 




n 


sinh 01 


cosh 01 


a 


n 


cosh ill 


"J*" J 

n ' 

■ 


where 


n* = (to*)^ - (k^ + H^)N^ 
n ' n ' n 


The boundary conditions (3.31) applied to Equation (3,41) yield 

+ Jq VP 


A = 
n 


k nV 
^ 1 _ 

J?(k ) ol*^ •^^'*' ctho)**' - 01 n cthoi 
0^ n n n n n n n 


kVH^ 

1 




Jt(k ) 0 )’*' a**" tghoi*^ - 03 a thhcii 

O n nn^’n nn n 


■ijSir%r-iiL:^TST 


7aro:z:x 




r:^E 


(3.40) 


(3.41) ‘ 


(3.42) 


(3.43) 


(3.44) 




1 n“l 


m'JI' cthti)' - ( 1 ) n cthw 
n n n. n n n 

siiih 5 
n 

sinh u ?\ 
”1 


*1- 

sltth ti} 
n 

sinh (d / 
n 



/ cosh Id*' ? ctzah u ? 

/ n n 


u n tghtii - tii S2 tghct) vcosh cc cosh u 


P) 


to fi tghcj ^ (0 n tghw 
nn^n un^n 


cosh to cosh to 

n n 


(3,48) 


The remaining dimensionless centrifuge fields E = ~ V4/E and J = J/J 

o c 

are given in terms of the solutions for $(p,C) and V(p,c): 


E = - 3$/9p , E- ■= 0 , E^ “ "N 3$/a? 
P 05 


<3.49) 


T ss ^ 

p i + 


/ t. TT\ 7 _ tOT ^ I *tA 7 

3p * “^0 TTw^ 3p 


where E s $ /R and J = o4 /R [Eq. (3.13)], 
000 o 00 


(3.50) 


If the cathode is in the plane z = ~c(? = -1) and the anode is in 
the plane z - +c(? = +1), then the reference fields and [Eq. (3.13)] 
are negative, since I < 0. The resu3.ts are also applicable to the case 
where the anode is in the plane z = -c(5 = ~1) and the cathode is in the 
plane z = +c(^ + +1). In the latter situation, the reference fields V 


and are positive, since I > 0. These explanations hold for magnetic 
fields pointing in tlie positive z-direction, >0; changes its 
sign if B < 0 [Eq. (3.13)]. 
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3.4. Numerical Illustrations and Results 

As an illustration the radial (p) dependence of the dimensionless 

centrifuge fields V(p,c), $Cp,C)» E Cp,C)» E (p,?), and J (p,0 has 

P 5 P 

been calculated for I < 0 in the cross-sectional planes Z ^ -0.99 

(cathode region), ? =* 0 (central region), and C = 40.99 (anode region) 

based on Equations (3. 47) -(3. 50) . The remaining fields J_(p,c) and 

J-(pj?) are proportional *to J (p,C) and E (p,s), respectively 
6 P 5 

[Eq. (3.50)]. The characteristic dimensionless magnetic interaction 
numbers are treated as parameters: 

MT *= 1, 10; H = 1, 10, 100, 

The geometry parameter N is taken to be N = 1 so that M = (1 + w^t^) 
corresponding to = c [Eq. (3.20)]. The radial positions of the 
cathode and anode are assumed to be 

p = 0.01 (R = 0.01 R ) ; P. = 0.9 (R. = 0.9 R ). 

“ — O 'r "TO 

The dimensional fields are negative everywhere where the dimensionless 
fields are positive, and vice-versa [Eq. (3.12)] since < 0 and 
< 0 for I < 0 [Eq. (3.13)]. 

i) Ce.rvtJtaZ Reg-ion, ? = 0: In the Figures 3. 2-3. 6, the aximuthal 

velocity field V(p,0), the electric potential $(p,0), the radial and 

axial electric fields E (p,0) and E (p,0)®J (p,0) , and the radial current 

P 5 5 

density Jp(p,0) are represented versus p, td.th (aiT,H) ® (1,1), (1,10), 
(1,100), (10,1), (10,10), (10,100) as parameters. It is seen that |v[ 
increases considerably at any point 0 < p < 1 if either H or ut are 
Increased. In the region p a 0 sufficiently close to the axis, j4»l, 

•j 

|E |, |e - 2n"'‘‘'|, and [j 1 increase with increasing H or oiT, The field 
P Q P 

distributions move towards the axis p * 0 as wx becomes larger. The 
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Fig. 3.5. E_(p,c) versus p for 5 = 0 , and (o)x,H) = (1,1) £0 (10,100) 
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Big. 3,6 



Jp(P»C) versus p for C =■ 0, and <(ijt,H) = (1,1) to (10,100) 
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"hump" developing at p “ 0.9 (Figs. 3. 4-3. 6) with increasing oit shows 
the influence of the ring anode (p * 0.9, 5 “ +1) in the plane ? * 0. 

ii) Ccuthode. Region, c ® ~0.99: The Figures 3.7-3.11 show 

V(p, -0,99), 4(p, -0.99), E^(p. -0.99), E„(p, -0.99)«J_(p, -0.99), and 

J (p, -0,99) versus p with (tuT,H) » (1,1), ...(10,100) as parameters, 

P 

These fields increase in intensity at any point 0 < p < 1 if H or wt 
is increased. Since the ring cathode is at p = 0.01 (c = -1), the 
field distributions are closer concentrated at the axis p = 0 than those 
in the plane ? = 0 (Figs. 3. 2-3 .6). Note that the plasma rotates only 
in the region p = 0.1 with a significant velocity, since the Lorentz 
force decreases rapidly with Increasing p 1. A comparison of 

the corresponding fields in Figures 3. 2-3. 6 and Figures 3.7-3.11 
indicates that the discharge spreads slightly in radial direction with 
increasing -1 < ? <_ 0. In particular, an Increasing radial section of 
the plasma rotates with a significant speed as -1 < ? ^ 0 increases. 

lii) Anodz Reg-tPft, 5 » +0.99: In the Figures 3.12-3.16, 

?(p, +0.99). $(p, *H}.99), E (p, 40.99), E (p. 40.99)«J (p, +0.99), and 

Jp(p, +0.99) are plotted versus p with (tuT,H) = (1,1),... (10, 100) as 

parameters. The dependence of these fields on H and mr is as in the 

previous cases for C “ 0 and 5 ** -0.99. The velocity distributions are 

fully developed nearly through the entire chamber across section 

0 < p 5 0.9, since the Lorentz force -J^B^is strongest in the vicinity 

p “ 0.9 of the ring anode, p “ 0.9(? ■ +1). As a result, a thin and 

steep boundary layer exists close to the cylinder wall (p *= 1) with 

backflows at sufficiently small tor-values (Fig, 3.12). The radial 

distributions of $, E , and J (Figs. 3.13-3.16) clearly indicate 

P ? 5 P 
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:Jv i, 


Ctiat« In the plane 5 = +0.99, the electrical discharge has shifted to 
the region p - 0.9 due to the influence of the nearby ring anode, 
p_j^ * 0 . 9(5 “ +1)* This shift occurs first slowly in the region 
-1 < 5 < +1 - A 5 , and then rapidly in a relatively thin layer i5 << 1 
close to the anode plane 5 = + 1 . 

It is remarkable that the discharge remains concentrated in a 

radial region close to the cylinder axis with little radial spreading 

of the current density J, except in a layer close to the ring 

electrode of large radius (R. » R ) in which the radial current 

*r ^ 

component dominates the axial current component J^. This spatial 
concentration of the discharge is the more pronounced the larger H and 
WT. The speed of plasma rotation V(p,5) increases with increasing 
magnetic induction by orders of magnitude over the reference speed 
as the Figures 3.2, 3.7 and 3.12 indicate which show V(p,5) for 
increasing wT and H. The theoretical electric field and current 
density distributions are in qualitative agreement with experiments 
(Schwenn 1970) . 

The graphs in Figures 3.2-3.16 are based on the Fourier-series 
solutions, in which the first 100 terms were considered and the eigen- 
values k^were calculated up to the 10th decimal point. An even larger 
number of terms in the Fourier series solutions has to be taken into 
account if one wishes to compute (approximately) the centrifuge fields 
ectremely close to the ring cathode (p_ = 0.01, 5 = —1) and ring anode 
(p_j_ = 0.9, 5 = + 1 ) where 3§(p,5)/35 changes discontinuously with p 
doe to the electrode boundary conditions. The Fourier-series represen- 
tation may be highly unreliable for precise numerical work in the 
vicinity of a discontinuity due to Gibbs phenomenon, i.e., the overshoot 
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at a discontinuity. The Gibbs phenomenon could be suppressed 
drastically by the use of Lanczos convergence factors, which also 
accelerate convergence and are used as a means of data smoothing. 

Ihe system analysis presented Indicates that extremely high 
speeds of plasma rotation are obtainable already at moderate discharge 
currents I and magnetic Inductions presuming the magnetic inter- 
action numbers are not small, H > 1, wt > 1. As an example, consider 
a plasma centrifuge with 

" |lj = 10^ amp , jB^l = 10® Tesla, 

2 —1 
cf = 10 nflio/m , R = c *= 10 m. 

o 

Hence, by Equation (3,13) 

V “ Ic/2TroB = (5 /tt) x 10^ m/sec, 

• o o o 

and, by Figure 3.2, 

0[V_] = 0[V V] = 10^ m/sec , for lax <= 10, H = 100. 

0 o 

Speeds of plasma rotation V., which are by orders of magnitude larger 

than 10^ ra/sec, can be produced if the order of magnitude of the 

parameters car and H is increased. The viscous forces reduce, however, 

the speed of plasma rotation always in the layers close to the walls 

(z = +c; r = R ), 
o 






31^ 
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CHAPTER IV 

EFFECT OF INDUCED MAGNETIC FIELD ON DIVERGING PLASMA CENTRIFUGE 
AT ARBITRARY MAGNETIC REYNOLDS NUMBERS 

This chapter deals with the boundary-value problem for the partial 

differential equations, which describe the (azimuthal) rotation velocity 

and induced magnetic fields in the diverging plasma centrifuge with ring 

electrodes of different radii and an external, axial magnetic field. 

The closed— form solutions of the Fourier— Bessel series are obtained 

based on the magnetogasdynamic approximation for dense isotope plasma 

with negligible Hall effect. The electric field, current density, and 

velocity distributions are discussed in terms of the Uar *n»==»nTi number H 

and the magnetic Reynolds number R. For small Hall-coefficients, 

WT « 1, the induced magnetic field does not affect the plasma rotation. 

3 

The rotating plasma with speeds as high as 10 m/sec is obtainable at 
typical conditions that can be realized in the practical application to 
the Isotope separation. 

4,1. Boundary-value Problem for Velocity and Induced Magnetic Field 

The plasma centrifuge model under consideration is the same as 
depicted schematically In Figure 3.1 of Chapter III. The plasma is 
sustained by a discharge current I, which enters the centrifuge chamber 
of radius R^ through a ring anode of radius R_^ in the anode plans 


z B *fc and leaves it through a ring cathode of radius in the cathode 

r R • r R . 

plane z - -c. Accordingly, 2ir / ° J (r,z) rdr -if ° 6(r-R.) dr = I 

j Q Z J Q ± 

for the axial current density J in any plane -c < z < +c. The external 

magnetic field is axial and homogeneous = {0, 0,B}, In view of the 

o o 

symmetry of the system with respect to the axis r = 0, the plasma flow 

field is azimuthal, ^ = {0, Vg(r,z), 0}, so that V • ^ = 0, i.e., the 

plasma brfiaves incompressible. For negligible Hall-effect (djx«l) , 

Jg »= 0 and V X B = 0, J^} in accordance with Maxwell* s equation 

for the magnetic induction B. Hence, = 0 and = 0 because of the 

homogeneous boundary conrllcions for B and B . However, B-(r,z) ^ 0, 

'1C Z w 

since 

B-(r,z) _ = p l/2irR , -c < z < +c , 

0' * 'r=R o o — — 

o 

[rB_(r,z)] = P I 6(r-R )/2TTr , 0 < r < R 

r9r Q *'"* z=±c o ± — — o 


Since the induced magnetic field (B„) is azimuthal, the induced 
electric field (E ) is due to the rotation (V.) of the plasma across 

X u 

the external magnetic field (B ). The pressure distribution P(r,z) in 

o 

the rotating plasma is deteraiined by the r- and z-components of the 
magnetogasdynamic equation of motion [Eq. (2.1)], 



where 


J B 

3r z e* 


9z 


J B„ 
r 0 


(4.1) 


1 a 

r 9F<^V = Vz* “ Vr 


(4.2) 


The current density, J(r,z), and pressure, P(r,z), fields are 

-y. 

readily determined from the magnetic field B = {0, B„, B } and the 

O 


[ 

I 


i .lil J-] r LJ. 1 



1 




51 


velocity field V = {0, Vg, 0}, whereas the electric field is given by 

y ~y~ 

Ohm's law, E » - V x b + J/cr. In accordance with the equation of motion. 
Maxwell's equations and Ohm's law (Section 2,1), the plasma in the 
centrifuge with homogeneous magnetic field (Fig. 3.1) is described 
by the boundary -value problem for azimuthal velocity Vg(r,z) and 
azimuthal induction B„(r,z) fields: 


32v 




[i (rV )1 + 

3r ^r 3r ^ a 2 uu 3z » 

aZ O 


(4.3) 


a 1 a 

^ (rB„) 1 + — =— = - y oB — ^ , 
3r r dr 0 * 2 o o 82 

az 


(4.4) 


where 


Vg(r,z)r^j^ =0, -c £z £4c , 


Vo(r,z) , 
6' z= ±c 


- 0 , 


0 < r < R , 
— — o * 


(4.5) 

(4.6) 


and 


B0(r,z)j.^R = v^I/ 2 vR^ , -c<z<+c. 


1 3 




(4.7) 

(4.8) 


Equations (4.3) and (4.4) are the azimuthal components of the 

equations of plasma motion and magnetic induction = -y oVx(v x B) ] . 

o 

The boundary conditions (4.5) and (4.6) consider that the plasma does 
not slip at the walls r == R^ and z = ±c. The boundary conditions (4.7) 
and (4.8) follow from Maxwell's equations for a total discharge current 
of |l[ amps flowing from the ring anode (r=R^) to the ring cathode (r=R ) 
of vanishing radial width if I < 0 (Fig. 3.1). 


4.2. Analytical Solutions In Terms of Fourler-Bessel Series 


For physical and mathematical reasons (Section 1.3), It is 
suitable to formulate the boundary-value problem for the coupled 
plasma fields V_(r,z) and B.(r,z) in dimensionless form by introducing 

O w 

the dimensionless independent and dependent variables, 


p = r/R 

, 0 IP 1 1 , 

(4.9) 

0 

C “ z/c. 

-1 1 s 1 +1 , 

(4.10) 

V(p,c) = 

Vg(r,z)/V^, B(p,c) = Bg(r,z)/B^ , 

(4.11) 


where the reference values V and B are defined as (B = external 
>- o o o 

li^uctlon) , 

V = l/2irRB ac, B = B . (4.12) 

o o o o o 

In the dimensionless formulation, the boundary-value problem for 
the azimuthal velocity, V(p,?), and azimuthal induction, B(p,C), fields 
assumes the form 


3 rl 3 / t,m «-2 a^V H 3B 
.j^[-^(pv)] +N ----- , 


3p P op ^~Z 


ac‘ N 


^ av 

2 3? » 


where 


-11 4 1+1 . 


V(p,?)^^^ = 0, 0 IP 11 


(4.13) 


(4.14) 


(4.15) 


f 


(4.16) 


B(p,C)p^l = R, -lie 1+1 , (4.17) 

[pB(p,c)]^_±j = R fiCp-p+)/p, 0 £p £1 , (4.18) 

with 

H H (a/p)^B R , N S c/R , R = y I/2 ttR B = y pV c^O. (4.19) 
oo o o oooo< 


The Hartmann number H, N, and the magnetic Reynolds number R characterize 
the ratio of Lorentz to viscous forces, the geometry of the centrifuge, 
and the intensity ratio of the induced and external magnetic fields, 
respectively. Equations (4.15), (4.16) and (4.17), (4.18) are the 
homogeneous and inhomogeneous boundary conditions for the fields V(p,4) 
and B(p,^), respectively. The linear statement. 


B(p,0 = R[p + '1'(P»?)I , (4.20) 

reduces the Equations (4.14), (4.17) and (4.18) for B(p,?) to equations 
with a homogeneous boundary condition (4.22) for '?(p,4) 


3p 




. -2 _ „-2 8V 

+ N - N ^ 

a?2 


( 4 . 21 ) 


where 


np,C)p^3 = 0, -ll?ltl , 

5(P“P+) 

> 


(4.22) 


^ 4 : tpnp.c)],^, 


2 


0 < p < 1 


(4.23) 



v^Cp.O = J,(V) f^(5) , 


(4.24) 


¥ (P,?) = JiCfc P) g„(?) , 

II i. n n 


(4.25) 


where the eigenvalues > 0 are determined by the homogeneous boundary 
conditions (4.15) and (4.22) as the real roots of the transcendental 
equation. 


Jl(k^) = 0 , n = 1,2,3,... 


(4.26) 


Thus, the general solution of the coupled Equations (4.13) and (4.21) 
obtains by linear superposition as the Fourier -Bess el series. 


V(p,?> = E , 

n=l 


(4.27) 


"ifCpsO = 2 ^iCk^P) Sj^C?) . 

n=l 


(4.28) 


g’ ' - g “ -f ’ 


n “n n 


(4.30) 


By elimination. Equations (4.29) and (4.30) are reduced to decoupled 


differential equations of fourth order. 


I 

i I i 


Introducing Bessel's function J. (k p) of first order, partial 

n 

solutions of the coupled inhomogeneous Equations (4.13) and (4.21) are 
in the form. 


! I . I 


Substitution of Equations (4.24) and (4.25) into Equations (4.13) and 
(4.21) yields ordinary coupled differential equations of second order 
for f^(?) and gj^(?) : 


f " - k^ f = g' 

n n n n 


(4.29) 


\\ 


.. J .’; : . l .„ . i ,.., ■ 1 . 


. 1:11 


f » f 

n 


“ (2k^ + f * 0 , 

u ' n n n _ 

- (2k^ + H^)N^g' * + k^ g =0 , 
' n n * 


(4.31) 


(4.32) 


f (C)-„., 

n ?=±1 


(4.33) 




(4.34) 


as boundary conditions, by Equations (4.16) and (4.23), respectively. 

In deriving Equation (4.34), the Fourier-Dini series representation of 
the Dirac function [Eq. (3.32)] has been used. In addition to Equations 
(4,33) and (4.34), f^CC) and have to satisfy also the coupled 

Equations (4.29) and (4.30). With the four real roots of Equations 
(4.31) and (4.32) [f . g “exp(u^C )]5 

li 11. n 


4» 

a- H w , = (1) 

In n 2n u 


*3n = "“n ’ “4n = "“n ’ 


(4.35) 


where 


u" = 2“^ N{(2k^ -F H^)±[(2k^ + - 4k^]^}^ , 

n ' n ' n n^ 


(4,36) 


the general solutions for fjj(C) and of Equations (4.31) and (4.32) 

can be written as 

. sinh w C . cosh tn c 

- A+ 1- + 1- 

n n • 4| * xt * 'J 

sxnh (i) , cosh w 

HU /A 


(4.37) 


sinh u ^ 


5 cosh 111 5 
S_ q. b“ S- 


sinh ( 1 ) 


° cosh u ^ 


n 


n , H 
cosh 01 

T1 


sinh 03 x; 

+ c” ^ + d" 

^ slnh 01 ” 

n 


cosh (I) C 

n_ 


cosh 0 ) 

n 


(4.38) 


+ 4 - 

Only four of the eight integration constants A^, D“ are independent. 
Substitution of Equations (4.37) and (4.38) into Equation (4.29) and 
Equation (4.30) yields 


A" [(6i“)^ - = -hV tghw~ D* , (4.39) 

n 'n n n °nn* 

B“ [(w“)^ - kV]/tii" = cthoj" C"' , (4.40) 

n^'n n n nn’ 


+ -+2 22 .+ + + 

[(to;) - k^K ]/to; = - tghto; b; , 

-1- 49 994 44 

D“ ECto”) - k N ]/to = - cthfcj“ A , 
n^ u n n nn’ 


(4.41) 

(4.42) 


respectively. The coefficient determinant of Equations (4.39) and (4.42) 
or Equations (4.40) and (4.41) vanishes (condition for existence of 
nontrivial solution) , 

A* 5 I(u~)^ - kV]^ - hV(io-)^ = 0 , (4.43) 

n n n 

in agreement with Equation (4.36). From the latter or Equation (4.43) 
one deduces the relatioas. 


E<<)^ - = ±HH . 

u n 12. 


(4.44) 


vhlch simplify the left-hand sides of Equations (4.39)-(4.42) , 



, ■ ji.. ‘ ' £- 


I L- 


Application of the boundary conditions (4.33) to Equation (4.37) 


shows that 


+A = A , 
n n * 


-B «= +B H B 
n n n 


(4 .45) 


Substitution of Equation (4.45) into Equations (4.37) and (4.38) gives 


( sinh (I) t sinh u ^ 
H ^ 

sinh ut sinh o) 


( cosh to C cosh 0) e;\ 

a ^ J ^ 

cosh tij cosh bi / 
n n 


(4.46) 


, /cosh u) t cosh ill C’ 

\sxnh w sxnh to 
n n 


, /sinh to C sinh to 
_T, if + S 

nMV T,” 

\cosh (I) cosh ill 
n n 


osh ill C\ 
— 

inh to / 
n 

inh to 

— ^)’ 

osh ill / 


(4.47) 


the latter under consideration of Equations (4.41), (4.42) and Equation 
(4,44). Application of the boundary conditions (4.34) to Equation 
(4.47) yields, upon elimination 


^ ^n (cthw^ + cthoi^) JQ(hj^) 


(4.48) 


“ (tgho)^ + tghto~) 


(4.49) 


By combining Equations (4.46)-(4.49) the solutions for f^^C?) and 
g^(5) in final form are 




5.8 


f (C)/NH = - 


+ -^0<V+> 


(cthoj + cthti) )k -nv 
n n n U n 


j;(k ) \ 


slnih (0 ^ 
n 


sinh (D ^ 


slnh 03 


n 


sinh 03 




n. 


J_(k p ) - J^(k p J 
0 n - 0 n*^+ 


(tgho)**" + tgho3 )k JftCk ) Vcosh o)"^ 
“n “n nO n n 


t( 


cosh 03*^5 
n 


cosh 03 ^ 


cosh 03 


f ) 


g^(0 = + 


J„(k p ) + JnC^t Pl) /cosh 0 ) 5 cosh oi Z 
0 n*^- 0 n + / n , n. 


+ - 2 \ 

(ctho3^ + ctho 3 ^)k^ J^(fc_) \sinh o) 


''A' 

II n U n 


sinh £i3 




Jq<V-> - Jq V+> / 

(tghio^ + tghai~)k J?(k ) V 


sinh 0 ) ? sinh to ? 
?^ + — " 


n 


n n O n 


cosh 0 ) 


cosh 0 ) 




(4.50) 


C4.51) 


Below, also the ^-derivative of Is required, which is given by 




el<?) *= + + _ 2 

(ctluo + ctho3 )k 


n 


n n u n 


( . Sinn (0 c Sinn w ?\ 

w 4- to J 

, ° sinh 03 ^ sinh 03** / 


n 


Jq<V-> " V . - “n^ \ ,, 

—y ( 03 V- + 03 ^ 1 . (4 . 

J„(k ) \ ”^ cosh 03 ^ cosh 03 / 


(tgho3^ + tgho3^)k^ Jq 
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n 


In terms of f (^) , g (^) , and g'(4)i the solutions for the 
n n n 

dimensionless fields V = {0,V,0}, B - {0,B,1}, J = {J ,0,J } and 

P 5 

^ *» {Ep,0,E^} of the plasma centrifuge are by Equations (4.20), (4.27) 
and (4.28) 


V(p,0 = 2 
n=l 


J (k p)f (O 
in n 


(4.53 


B(p,C) = R[p + Z J,(k o)g„(?)l . 


(4.54 


l.'i . 1 . i : .X. 1 ■ I 1 :^ 1 ^ 
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Jp<p,c) = - N"^ z J (k^p)g»(?) 

n«l 


J-Cp#5) = 2 + Z k J (k p)g (c) 
? n*=l ° ° 


(4.55) 

(4.56) 


Ep(p,?) ** - V(p,c) + N Jp(p»?), E^(p,c) « NJ^(p,5). (4.57) 

The reference values and for V(p,^) and B(p,c) are defined in 

Equation '4.12). The dimensionless fields J (p,C) and E (p,C) are 

P P 5 S 

normalized with respect to 

J = 1/2ttR^ , E H V B = l/2irR c . (4.58) 

0 0*000 0 


It can be proved that the solution V(pa?) derived here from the 

boundary-lvalue problem with a significant induced magnetic field B_ 

'remains valid even for the boundary-value problem in Chapter III as 

long as the Hall effect is negligible. Let V^j.j.(p,?) and V^^Cp,?) 

designate the solutions for the dimensionless azimuthal velocity in 

Chapter III for ut « 1 and the present chapter, respectively. The 

Ohm’s law for mx « 1, J = - + V_B , and Faraday's law, 

’ r 3r e o •' 

- T — = poj , give the following relation in dimensionless form 


3p 


+ V. 


Ill 


R 3? 


(4.59) 


Substitution of this relation into Equations (3.14), (4.13) and boundary 
conditions for and yields a new boundary-value problem for 

W(p,0 5 Vj.jj.(p,C) - Vj^(p.C): 


3 rl 3 / ttm J. «“2 3^W „ 

^!-s^(PW)l +N -I =0 • 


(4.60) 
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vbere 

W(p,5)p^l “ ° 

W(p,c)^.+3^ = 0 , 0<p<l . (4.62) 

The general solution of Equation (4.60) has the form 


W(p,C) = 2 J,(k p) [A^ coshCNk C) +B slnh(Nk ?)] , (4.63) 

n n n 


where the integration constants and are determined by boundary 


conditions (4.62): 


A cosh(Hk ) + B sinh(Wk ) “ 0 , 
n n IX IX 


(4.64) 


A cosh(Nk ) - B sinh(Nk ) = 0 . 
c n n n 


(4.65) 


However, W(p,c) should have only trivial solution (W = ~ 

since the coefficient determinant for A^ and of Equations (4.64) and 
(4.65) would not vanish, i.e.. 


cosh(Nk ) 
n 

cosh(Nk ) 
n 


sinh(Nk^) 

- sihh(Kk ) 
n 


^ 0 


( 4 . 66 ) 


This means that V^jj(p,?) [Eq. (3.47)] for tor = 0 should be identical 
to Vjy(p,?) [Eq. (4.53)]. 
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4«3. Numerical Illustrations and Results 

As an illustration, the radial (p) dependence of the dimensionless 
discharge fields VCp.C), B(p,?), E (p,5), J (p,C)» and J (p,?) has been 

MM ^ 

computed for I < 0 in the cross-sectional planes ? = -0,99 (cathode 
region), ? » 0 (central region), and ^ = +0.99 (anode region) based on 
Equations (4.53)-(4.57) . The remaining field E (p,c) is proportional 
to J^(p,^) [Eq. (4.57)]. The characteristic dimensionless magnetic 
interaction number H is treated as a parameter: H = 1> 10, 100. The 
geometry parameter N = taken to be N - 1 corresponding to 

^ c [Eq* (4.19)]. The radial positions of the cathode and anode are 
assumed to be 


p “ 0.01 (R = OcOl R ) ; p. = 0.9 (R, - 0.9 R ) . 


With the exception of — B B, the dimensional fields are negative 

9 o 

everywhere where the dimensionless fields are positive, and vice-versa 

since V < 0, J <0 and E < 0 for I < 0 [Eqs. (4.12), (4.58)]. Note 
o o o 

that the ma^etic Reynolds R in Equation (4.19) is defined to change 
its sign with the sign of „ 

The Equations (4.53)-(4.57) indicate that the velocity field 

V(p,C)» the current density field J _(p,?), and the electric field 

P» 4 

E (p»?) are independent of the magnetic Reynolds number R, whereas 
P»4 

the induced magnetic field B(p,?) is proportional to R. This is due to 
the azimuthal direction of the induced magnetic field B(p,?), which is 
parallel to the velocity field V(p,5) of rotation. Accordingly, the 
plasma fields V(p,C), B(p,C)/R, J _(p,C), and E (p,^) depend only on 
the Hartmann number H, presuming that the Hall effect is negligible 
(wT « 1) . 
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i) Ce.ntfLOZ Refl/COn, 5 = 0 : In Figures 4. 1-4 .5, V(p,0), 

[B(p,0) - Rp]/R, EpCPjO), Jp(p,0), and J^(p,0) « E^(p,0) are shown 

versus 0 <_ p £ 1 with H = 1, 10, 100 as a parameter. It is seen that 

jvj increases considerably at any point 0 < p < 1 as H is Increased. 

Similarly, (B - Rp)/R and the sources J ^ of the magnetic induction 

increase in intensity within the main central region 0 < p < 1 - ip as 

H is increased. For large values H > 10, B and J decrease in the 

~ P » ? 

wall region Ap = Ap(H), so that the electrical discharge becomes more 

concentrated in the center 0 < p < 1 - Ap of the centrifuge. The 

intensity of E increases uniformly in the region 0 < p < 1 as H is 
P 

increased, while E “ J . 

il) Cathode. Regtonf c ~ -0.99; The Figures 4.6-4.10 show 
V(p,-0.99), [B(p,-0.99) - Rp]/R, E (p,-0.99), J (p,-0.99), and 

H r" 

J^(p»“0-99) “ E^(p,-0.99) versus 0 ^ p ^ 1 for H = 1, 10, 100. The 
fields V, E , and J increase in intensity at any point 0 < p < 1 

P* b P J b 

with increasing H, whereas B/R decreases in 0 < p < 1 with increasing H. 

Since the ring cathode is at p_ = 0.01 (? = -1), the field distribu- 
tions are more closely concentrated at the axis p = 0 than those in the 
plane ? = 0 (Figs. 4. 1-4. 5). Note that the plasma rotates only in the 
region p = 0.1 with a significant velocity, since the Eorents force 
decreases rapidly with increasing p 1. 

iii) Anode. Reg-ton, 5 = +0.99; The Figures 4.11-4.15 present 

V(p,-H).99), [B(p,+0.99> - Rp]/R, Ep(p,+0.99), Jp(p,+0.99), J^(p,+0-99) 

« E^(p,+0.99) versus 0 <_ p <_ 1 for H = 1, 10, 100. The velocity field 

is fully developed nearly through the entire centrifuge across section 

0 < p < 0.9, since the Lorentz force -J B is strongest in the vicinity 
- p o 



Fig. 4.1. V(p,0 versus p for ? = 0, and H = 1, 10, 100 
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Fig. 4.10. J Cp,?)- 2 versus p for 5 = -0.99, and H = 1, 10, 100 












Fig. 4.15. J^(ps?)“2 versus p for ? = +0.99j and H = 1, 10, 100 
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p “ 0.9 of the ring anode = 0.9(5 ® +1). As a result, a thin and 
steep boundary layer exist? close to the cylinder wall (p * 1) with 
plasma counter-rotation at sufficiently small H-values. The radial 
distributions of B, Ep ^ clearly indicate that, in the plane 

5 “ +0.99, the electrical discharge has shifted to the region p = 0.9 
due to the influence of the (nearby) ring anode at p^ ** 0.9(5 +1) • 

In the graphical illustrations, the cathode radius R_ was chosen 
to be small compared to the anode radius R, to ensure a large angle 

"T 

between the current field lines J(r) and the external magnetic field 

i.e., a significant Lorentz force. A comparison of the Figures 

4.1 and 4.6 with Figure 4.11 indicates that this choice of electrode 

radii results in a radial boundary layer of large width and low velocity 

in the lower half -c ^ z < 0 of the centrifuge. Hence, R « R 

(or R_ » R_j_) is not the best choice for a centrifuge of maximum 

efficiency.' Figure 4.11 demonstrates that a velocity rofile rising 

uniformly with radius r and decreasing rapidly in a steep boundary 

layer of narrow width 4r, is obtained by using a cathode and an anode 

of the same radius R_ = R_^^ < R^, which is nearly as large as the 

centrifuge radius R^. Although R ~ R. in this case, the current field 

lines J(r) intersect with B at a sufficiently large angle 3f(^B ) 5 ^ 0 

o o 

due to the repulsion of the current filaments. As a result, a net 

Lorentz torque restilts for a centrifuge with R = R . which is still of 

the same order of magnitude as for a centrifuge with R « R, (presuming 

that I, and B , c, and R are the same) . 
o o 

The centrifuge analysis presented indicates that extremely high 
speeds of plasma rotation are obtainable as shown in Chapter III at 
moderate discharge currents I and magnetic inductions B^, presuming the 


Hartmann number H is not small, H > 1. As an example, consider an 
isotope centrifuge discharge with the same values of |l|, o, 

and c as those in Chapter III. 

Hence, by Equation (4.12) 

V ® I/2trR B crc “ (5 /tt) x 10^ m/sec, 

and, by Figure 4.1 

0[Vg] = 0[V^V] = 10^ m/sec, for H = 100. 

Based on these examples, one can assume with some confidence that 
hl^"Power plasma centrifuges are technically realizable employing dense, 
collision-dominated isotope plasmas. The proposed high-density plasma 
centrifuge would use arc plasmas at pressures of about one atmosphere 
so that the isotope masses separated are increased by orders of 

h 

magnitude. The large Hartmann numbers H - (cf/y)^ for high 

speeds of isotope rotation are achievable because of the (relative) 

small viscosity y and large conductivity c of gaseous plasmas. Speeds 

of plasma rotation, which are by an or der-of -magnitude larger than thosr 

in the above examples, can be achieved at realistic Hartmann numbers H, 

11 

Since m = 1.76 x 10 B sec., the Hall-effect is insignificant in dense 

-12 

plasmas for B = 1 Tesla as long as T < 10 sec. In general, the 
.Hall-effect increases the speed of plasma rotation for 6 jt > 1, i.e., in 
plasmas of lower density as shown in Chapter III. In developing a 
plasma centrifuge, therefore, apparently a trade-off between isotope 
density and rotation velocity has to be made. 
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CaiAPTER V 

PLASMA COUNTER-ROTATION IN MULTI-DISCHARGE CENTRIFUGE 

This chapter is concerned with a plasma centrifuge between two 
ring electrodes embedded in the mantle of a cylindrical chamber, in 
which the plasma in the anode and cathode regions rotates in opposite 
directions under the influence of a spatially converging and diverging 
current density and an external axial magnetic field. The associated 
boundary-value problem for the coupled partial differential equations 
describing the azimuthal velocity and radial current density fields is 
solved in closed form. The difficulties associated with the complex, 
inhomogeneous boundary conditions are overcome by means of Fourier 
expansions for Bessel functions of complex argument. The velocity, 
current density, induced magnetic induction, and electric fields are 
presented for typical Hartmann numbers, magnetic Reynolds numbers, and 
geometry parameters. The discharge is shovm to produce anodic and 
cathodic plasma sections rotating at speeds of the order 10 m/sec for 
conventional magnetic field intensities. Possible application of the 
magnetoactive discharge as a multi-discharge plasma centrifuge for 
Isotope separation is discussed. 

5.1. Model for Multi-discharge Centrifuge 

The plasma centrifuge to be studied herein exhibits the interest- 
ing effect of plasma counter-rotation, i.e., the plasma in the anodic 
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and cathodic half-spaces rotates In opposite directions . As depj ' ^ed 

schematically In Figure 5.1, the centrifuge system consists of an 

electrically Insulating cylindrical chamber of radius with end walls 

at z>=±L. A perfectly conducting ring anode (r=*R^, z ==-c) and ring 

cathode z=+c) are embedded in the cylinder R^ (eventually in 

form of thin, "hollow” slit electrodes) . The electrodes are placed far 

from the end walls (c « L) to reduce velocity losses due to the end 

plates. The plasma is produced in the space -c < z<+c, 

through a gaseous discharge resulting in a curved current density 

distribution J(r,z) which intersects the axial, homogeneous magnetic 

field B applied from outside. The J x B force rotates the plasma 
o o 

counter-clockwise in the s jode region -c < z < 0 and clockwise in the 
cathode region 0 < z < +c (Fig. 5.1), since the J-lines converge for 
z < 0 and diverge for z > 0. In the central plane z^, the plasma is 
at rest. 

The purpose of the investigation is to evaluate theoretically the 
electromagnetogasdynamlcs of plasma counter-rotation as a contribution 
to the physics of plasma centrifuges. Furthermore, it appears that 
this type of centrifuge might be useful for isotope separation, in 
which the heavy isotope would be enriched in the anodic and cathodic 
'‘layers" z - ±c and the light isotope would be enriched in the central 
layer z « 0 . The end spaces -L ^ z < -c and +c < z ^ +L would serve as 
reservoirs for the unseparated isotope mixture. For a proper choice of 
the geometry parameters, R^, c, and L, well-developed azimuthal velocity 
distributions can be expected. Iijf view of the circumferential electrode 
arrangement, a large number of such centrifuges can be set up in a long 
insulating cylinder to rotate large volumes of isotope mixtures, as 




Fig. 5.2. Scbeme of multi-discharge centrifuge. 

5.2. Boundary*>value Problem for Velocity and Radial Current Density 

The steady-state rotation of the plasma centrifuge shown in 
Figure 5.1 is theoretically investigated based on the magnetogasdynamlc 
equations (Section 2.1) for dense plasmas. Laminar flow is assumed and 
(conceivable) secondary flows superimposed on the main rotational flow 
are disregarded. Experiments Indicate secondary flows in the motion of 
incompressible fluids between rotating cylinders (Chandrasekhar 1961) , 
but secondary flows have not been observed in plasmas which rotate 
under the influencj^ of electromagnetic forces. In view of the symmetry 
of the centrifuge configuration with respect to z-axis, the plasma flow 
field is then azimuthal, V = {0, V (r,z), so that the plasma behaves 
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Incompressible (7*V=0) . It is assumed that the gyration frequency w 

of the electrons is much smaller than the collision frequency t ^ 

between electrons and plasma particles (o>t« 1) . In this case, the 

current density is of the form J « {J (r,z), 0, J (r,s)}, and the 

Hall-effect is negligible (dense plasmas of low ionization degrees) . 

The magnetic induction is of the form B = {0, B-Cr,z), B } in 

y o 

accordance with l-laxwell's equations and the homogeneous boundary 

conditions for B and B . 

r z 

The counter-rotating plasma centrifuge is described i»y the 
boundary-value problem for tSe azimuthal velocity Vg(r,z) and radial 
current density J^(r,z) fieldi:;: 


3^V B 

-SjA — (rV ) 1 + ^ ^ J 

^r -2 p r * 

dZ 


(5.1) 


a 1 a 

^^r air^^'^j^ ^ *** 0,2 *■ .,2 * 


(5.2) 


where 


Vg(r,z)^^^ =0, -L<2<+1, , 


(5.3) 


Vg(r.z)^^,j^ - 0, 0<r<R^ . 


(5 .A) 


Jj.(r,z)^^^j^ ” 2iiR“ ~ fi(z+c)l, -L £ z £ +L , 


(5.5) 




Equations (5.1) and (5.2) are the azimuthal component of the 
equation of plasma motion and the Induction equation combined with 
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V X f ** p J and 7«J = 0, respectively. The boundary conditions (5.3) 
o 

and (5.4) specify that the plasma does not slip at the chamber walls 

r “ R and z = ±L. The boundary conditions (5.5) imply that a total 
0 

discharge current I flows from the ring anode (z = -c) to the ring 

cathode (z = 4c) of vanishing axial width, Az -»• 0. The boundary 

conditions (5.6) consider that no radial current flows at the end 

plates at z = ±L according to Ohm's law, = o(E^4VgB^), since 

Vg(r,z)^^^j^ = 0 and \(r,z)^^^j^ = 0 by n x [f] =0. 

The remaining centrifuge fields are consecutively determined by 

using the solutions for Vg(r,z) and J^(r,z). The axial current density 

J (r,z) is obtained by the conservation equation for the electric 
z 

charge density (V • ^ - 0) : 

J^(r,z) = -J (rjpidz + c^(r) , (5.7) 

where the integration constant c^(r) is determined by the boundary 
conditions (no axial current flows into end walls) , 


J 

z 




= 0 , 


0 < r < R . 
— — o 


(5.8) 


The induced magnetic field Bg(r,z) is obtained from the z-component of 
Maxwell’s equation, V x B = u 


Bg(r,z) = C 2 (z)], 


(5.9) 


where the integration constant C 2 (z) is determined by the boundary 
conditions [equivalent to the boundary conditions (5.5) nince 
“3Bg/3z = yg*7j.» H(c-lzl) “ Heaviside step function]. 




p I 
o 


2ttR 


H(c-lzj), 


-L < z < 4-L 


(5.10) 


The electric field E(r,z) is given by Ohm’s law: 


I 

i 

I 

i 

i 

I 

i 

I 

i 


1 


I 


E(r.z) = 

Finally, the pressure distribution can be calculated from the r— and 
z-component of the equation of motion [Eq. (2.1)1 the plastma if 
desired . 


5.3. Fourier Series Solutions in Terms of Bessel fractions of 
Complex Arguments 


In order to solve analytically the boundary- value problem for the 
coupled plasma fields Vg(r,z) and J^(r,z), it is convenient to 
formulate Equations (5.1) -(5. 6) in dimensionless form by introducing 
dimensionless independent and dependent variables. 


and 


p - r/R^, 0 p 5 1 , 

5 ® z/c , ± 


L = L/c , 


( 5 . 12 ) 

( 5 . 13 ) 


V(p.?) « Vg(r,z)/V^, dp(p,0 = J^(r,z)/J^ , 

where the reference values V and J are defined as 

o o 

V = I/2irR cBa, J =ovB = l/2vR c 
o oo o oo o 

Thus, rhe boundary-value problem defined in Equations (5,i)-(5,6) 

becomes for V(p,c) and Jp(p,C): 


op P 3p g^2 P 


3p P op P a-2 




( 5 . 14 ) 


( 5 . 15 ) 


( 5 . 16 ) 


( 5 . 17 ) 


where 


V(p,C)p^j^ *0, -L < ^ <+L , (5.18) 

V(p,?)^^+2. " 0 1 P 11 » C5.19) 

Jp(p,C)p^l = -5(5+1) + 5(5-1), -L < 5 < +L , (5.20) 

Jp(p,5)^=+f^ = 0, 0 <^p _< 1. (3.21) 

The dimensionless parameters N and H are defined by 

N E c/R , H = (a/u)\ R . (5.22) 

o o o 


In view of symmetric geometry of the centrifuge configuration and 
asymmetric boundary conditions (5.20) and (5.21), J^(p, 5 ) is asymmetric 
about 0, i.e., Jp(p,5) = -Jp(p,“5)» Accordingly, the radial current 
densities at (r, ±z) have opposite directions, but same magnitudes. 
Consequently, as a result of the radial current densities normal to an 
axially applied magnetic field, the Lorentz forces set the plasma in 
rotation with an azimuthal velocity field which is also asymmetric, 
i.e., V(p,c) = -V(p,-0. 

In accordance with the above conclusions and the boundary 
conditions (5.18)-(5.21) , the general solutions of the coupled partial 
differential equations (5.16) and (5.17) are sought in the form of the 
Fourier series: 

00 

V(p,5) = S f (p)sinX^C, (5.23) 

n=l "" “ 

CO 

J (P,5) = s g„(p)sinx c. 


(5.24) 


where 


(5.25) 


\ = nir/i. 

n 

The p-dependent functions f^^Cp) deteriulned by the follow- 

ing coupled ordinary differential equations and boundary conditions: 


d rl d , , M - H^a 

d? '7rp- ‘I’V’ - 'T> “ ®n’ 


(5.26) 


ri Cl , V ■■ / n.2 _ / 2. - 

[-^ (Pg^)l - (r^-) g„ - (;,^) V 


rl _d 

dp dp 


(5.27) 

acj-j-co jLii t-iie xnrervax ~L ^ < +L , 


-6(5+1) + 6(5-1) =» (2/t) Z sljik sinX 5 . 

__T n n 
n— i 


(5.30) 


By eli m i n ation. Equations (5.26) and (5,27) are reduced to decoupled 
differential equations of fourth order. 


dp' 

rd rl d . 


:^)2>2 

N'' ^ 

^2,2 


2 * 

(5.31) 

3 

^ ®n * 

(5.32) 


The general solucions for fj^(p) ®®d Sjj(p) of the identical biharmonic 
equations (5.31) and (5.32) are 


e__Cp) - 

where A^, and arc integration constants. d 3 ^(**^jjp) ®®d d 3 _(*'*P) 

are complex conjugate Bessel functions of the first kind of order 1, 
since and k* are complex conjugate eigenvalues given by 


n 


K* 


n 


r^n, „ +i4> 

r^n, ^n -tim% 


(5.35) 
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where 


« N 


li 

% arccosi s-^^ 5 — 1 = h arcsin[ ^ 2 J<^ 

ncx^/hW)*^ cx;/nW)^ 


(5.36) 


Furthermore, the solutions (5.33) and (5.34) have also to satisfy 
original coupled equations (5.26) and (5,27) and boundary conditions 
(5.28) and (5.29), i.e., the integration constants are interrelated by 


sinX . 

. _ ^ MH n_ _ 

n Hx J,(k ) n * 
n 1 n 


sinX 

f* — S = rj* 

n LJ- (ic ) n * 
1 n 


(5.37) 


By combiitlng Equations (5.23) and (5.24), Equations (5,33), (5.34) and 
(5.37) and noticing that solution for V(p,5) and 

obtained in the final form 

ea 

V(p,?) = -2NHL”^ S X sinX^sinX^? Im[Jj^(K:^p) /J^( k^) ] , (5.38) 

n=l 


J (p»5) = S sinX^sinX^?Ra[J^(K^p)/J^(tc^)] , 

P n=l 


(5.39) 


where Re[Jj^(Kj^p)/Jj^(K:^) ] and Im[J^(ic^p)/Jj^(K^^)] refer to real and 
imaginary parts of Jj^(ic^p)/J^(ic^), respectively. 

The dimensionless expressions of Equations (5,7)“(5.11) together 
with the solutions (5.38) and (5.39) for V(p,?) and Jp(p,?) yield the 
remaining dimensionless discharge fields J^(p,?), B(p,C) and ^(p,?) 
which are normalized with respect to ^ ^^o^o ^o “ ^o^o* 
respectively: 


■ t «■ 1"' 


90" ; 


CO 

J (p,c) = 2HL"^ Z X“^sin\^[cosX^C-(-l)^] Re[K^jQ<(c^p)/Jj^(K^) ] , (5.40) 


n=l 


CO 

B(p,C) = 2RL"^ 2 x"^sinX^tcosX RaU^ (k p)/J, (k J1 , (5.41) 

-XI Xi n xn**-ii 

n»j. 


Ep(p,?) = Jp(p,C) “ V(p,5) , 


(5.42) 


E^(p,C) = J^(P,5) . 


(5.43) 


3n(.K p) is a Bessel function, of the first kind of order 0 with complex 
0 n 

argument, and R is the magnetic Reynolds number. 


R = y oV c = u I/2 ttR B . 
o o o o o 


(5.44) 


From all plasma fields, only the induced magnetic field B(p,0 depends 
on the magnetic Reynolds number. It is also noticed that Equation (5.40) 
satisfies the integral condition. 


J (p.Opdp = NHd-UI) , 


(5.45) 


vhich is proved with the help of the following Fourier series 
expansions in tb'i interval, -L ^ ^ 

5 » “2 £ (-1)^ X"-^sinX^5 , 

- n XI 

n=l 


(5.46) 


H(l-Ul) « L“^(l + 22 X“’^sinX cosX c) . 

^ —-ill o. n 

IL^X 

Equation (5.45) is rewritten in dimensional form as a non-vanishing 




integral for any cross section -c < z < +c, 

.R 

2w / ° J (r,z)rdr =* lH(c- |zj) , 

•'n Z 


(5.47) 


which shows that the plasma is sustained by the total discharge current I. 









Fig. 5.4* V(p,5) versus p for ? 







ii) CuAJfOVt dzn&^y J (p,?) [figs* 5.5-5.81; The 

P j ^ 

electric discherge becomes more closely concentrated near the cylinder 


wall as H is increased. Backward current density < 0) appears to 
be abundant near the cylinder wall for [ C [ > 1 as either H or N is 
Increased . 

ill) Induced magnctLc ^Zctd B(p,c) [Figs. 5.9 and 5.10]: The 

Induced magnetic field is denser near the cylinder wall at p = 1 than 
the axis at p = 0 as H is increased. B(p,C) widens over Z as either H 
is increased or N is decreased. For large N = 10 and small H < 10, B 
grows almost linearly with p in the interelectrode space, -1 < ? < 1, 

iv) ELejCtnXc (^leMh E ^ [Figs. 5.7, 5.8, 5.11 and 5.12]: The 

P> ^ 

intensity of E grows rapidly near the ring electrodes (p = 1, ? = ±1) 

and the r-iaximum of E shifts to the cylinder wall as H is increased. 

P>? 

E extends more widely along the ^-axis as either H is increased or 
N is decreased. 


The graphs in Figures 5.3-5.12 are based on the Fourier series 

solutions for the Bessel functions of complex argiaments, which are 

summed up to n = n terms such that n satisfies 
max max 


If(p»?)^ -5 ” 1 — lO”^, The numerical values 

max max max 

of complex Bessel functions Jq and in the series solutions are 
computed based on the algorithms by Gautschi (1964) . The Gibbs 
phenomena at discontinuities and the tendency of oscillation are 
suppressed by using the Lanczos convergence factors (Arfken 1970), 


3in[mr/(n^^ + 1 )] 
"’'/("max 


2 j m * * t 


max 


(5.48) 
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Fig. 5.8. J^(p,c) versus p for C “ 1.1, 

H « 1, 10, 100. 
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Fis. 5. 10, B(p, 5)/R versus p for C 1.1. 0,9, 0, N « 1, 10 and 
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The grephical presentations of the plasma fields indicate that 
it is desirable for the ring electrodes to he located sufficiently 
apart in distance compared to the cylinder radius (N “ ^ 1) » to 

ensure a significant Lorentn force and plasma rotation. As long as the 
end vails of the discharge chamber are placed sufficiently far from the 
ring electrodes (L » L/c » 1), velocity losses due to the boundary 
layers at the end walls are insignificant, llie proposed centrifuge 
scheme results in supersonic rotational plasma velocities (which are 
not affected by the induced magnetic field) for moderate flow numbers 
H, N» L that are realiaable in practical applications. For example 
|l| « 10^ amp, » 10® Tesla, o » 10^ mho/m and c « lOR^ ■ lO"^ m, 
the speed of plasma rotation is 10^ m/sec in order of magnitude by 
Equation (3.15) and Figures 5,3 and 5. A for U » 100 and U » 10, 

In a practical centrifuge design for isotope separation, the 
multi-discharge centrifuge would be located at some stage in a cascade. 
Hxe enriched ai\d depleted isotope streams are introduced at one end of 
each centrifuge stage and removed at the other end. As long as the 
inflow of the isotope mixture and the reiwoval of separated ions occur 
at a sufficiently slow rate, these flows can be neglected in the 
analysis of the plasma rotations. 
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APPLICATION TO ISOTOPE SEPARATION 

As shown in the previovis chapfcei?s» the plasma centrifuge using 
electromagnetic forces permits to generate speeds of rotation of the 
order 10^ m/sec for ordinary steai, -state conditions, which are by two 
orders of magnitude larger than the achievable speeds of about 
400 m/sec in the mechanical centrifuges. As a result of the high speed 
of plasma rotation, the plasma centrifuges could be xxsed with 
advantage for the separation of isotopes. The spatial separation of 
the isotope ion and atom components according to their particle masses 
is ensured by the strong centrifugal force aotiug on the isotope 
elements at those high speeds of rotation. Two different kinds of 
isotope components got different rotational velocities according to tho 
mass due to tho intorcomponent friction force. This frictional force 
together with the magnetic field results in a radial motion of tho 
isotope components, tho velocity of which is negligible compared with 
the velocity of rotation. The heavier components move outwards and the 
lighter inwards in the plasma centrifuge, Tlte diffusion for such a 
multicomponent plasma transverse and parallel to a magnetic field in 
the presence of a centrifugal force has been discussed by Bonnevier 
( 1966 )% For tho plasma centrifuges, it is possible to apply his theory 






, i 


■ I 




xm 

to get approxitt\atQ numerical results where the separation effect is 
taken into account. 

As an illustration to the plnsuux centrifuge, the isotope 
separation ratio et is calculated for two isotopes of lighter mass 
and heavier mass with tlie same charges, Tlie separation ratio at 
distances 0 <r;<R^ - <S, whore o is the viscous boundary layer 
thickness, is proportional to the power CWj-m^)Vg(r)/kT^ in both the 
ttcchonical and plasma centrifuges (T^' is the temperature of the 
isotope ions) , The separation ratio for the mechanical centrifuge 
is expressed by Cohen (1951) 


n.(r)/n.(r) <,m^-i\\.)V?(r) 

^ nTMTuTW “ 


] . 


( 6 . 1 ) 


while the separation ratio for the plasma centrifuge is given by 
Bonnevier (1966) 


n.(r)/n^(r) 

”p ' n. rO) 7 n.U>) * 


(m “mj)V*(r) R 

kT ' '' ^ r ^ ‘ 

o 


( 6 , 2 ) 


As a iJpeciflc examplo, loiuiKUn’ vm uranium plasma oentrifui^e 
containing the boavie* \sotop^' onu H i >ut the UghteriO ir 


to obtain eurieheJ I'' 


IS 


ol about <1 l voi’i tiu' nutural ubundimi't' et 


as nuclear fuel for reacioia. Kor tlte eiu ivhment ot nauuul utauium 
235 

with 0 , the centrifuge dischargx> \anild bo burnt either betwen sv-il id 

uranium electrodes or in an UFg atiuos^itere. Such uranium arc 
discharges can be operated at temperatures as low as ' 4500*K. In this 


r2? 


-20 


case, one has m23S''*"‘^3S " ^ “ 6.213 10 Joule. 

Hence, the isotope separation ratio for r =* R^/2 is 


tv" - V ' 
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-2J-. 


_;i, 1;-.-. 


) \ f: 


oi„ « 1,057 X 10® 
P 

for 


« 4.013 N 10® 

for 

VjCrt 

« 2,593 « 10^ 

for 

V(,(r) 


* 10 w/scic , 

•» 5 10^ ®/sQC! » (6.3) 

• 10^ ra/sec , 


It is seen that the plasma centrifuges could produce considerably 
larger separation ratios than of the mechanical centrifuges for 
the isjajtimum achievable velocity of about 400 m/sec. Consequently, the 
plasma centrifuge requires significantly fovfor stages in a cascade to 
get desired concentrations of the isotopes iii eoanparison with the 
mechanical centrifuges. 






I 



f 

-.,.4 


I 




107 


CHAPTER VII 

COMPRESSIBLE PLASMA CENTRIFUGE WITH SECONDARY FLOWS 

In this section, the theoretical analysis of the steady-state 

dynamics of a plasma centrifuge employing concentric cylinder 

electrodes and an axial magnetic field is proposed. The plasma is 

produced by a radial discharge of current density ? in the isotope 

mixture, and the rotation of the plasma is caused by the Lorentz- 

forces J X B , Based on the compressible magnetogasdynamic equations, 
o 

a mathematical method is proposed which permits to calculate the 
plasma fields, such as the velocity, mass and current densities, and 
the electromagnetic field in the centrifuge, as a superposition of 
primary and secondary fields. In this approach, the critical Reynolds 
number for the onset of secondary flows is determined as the eigen- 
value of the boundary-value problem for the secondary fields. Attention 
is given to the evaluation of the feasibility of plasma centrifuges as 
it is affected by secondary flows and viscous boundary layers. 

7,1, Plasma Centrifuge Model with Secondary Flows 

The subject of the consideration is a t\ro-dimensional theory for a 
plasma centrifuge with electromagnetic and viscous forces, the Hall 
effect and secondary flows. As a centrifuge model, the previously 
considered type is chosen with concentric cylinder electrodes and axial 
magnetic field (Fig. 7.1). This centrifuge geometry is symmetric with 
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segpQot to tUo mid tUo pldno e 0^ It should » tUovo£oi*o> mako 

tUo Ineluslon of soeot arjr flows Into tho analysis possible • Hiu 

isotope iid^ture is eontalned in the intereleottode space < r < 

iu £om of a plasma sustained by a tadial gas discharge between the 

cathode at i* « kg and the anode t « k^« Tlie tat’ial component of the 

ourrent density 5 fonts with the magnetic field a volume force 

^ the asiwuthal component (J ^ which produces the 

rotation in the eircuinferentlal direction “0* In the steady state, 

the J s U forces are balanced by viscous forces and Inertia forces of 
o 

the plasma motion. The resulting velocity field ^(r,a) of the plasma 
can be represented as a auperpositlmt of a primary rotation field (o) 
and secondary velocity fields In the r, G, and a -directions, 

^(r,a) ^ {0,v^g(r),0} *1* (v^Cr.a), VgCr,®), Vjj(r,a)} . 

Generally, the aaial extension of the centrifuge Is large compared 
to the radial extension Itg ~ 
treated as one-dimcnaional, e,g,, particular, 

this assumption is rigorous for a hypothetical centrifuge of infinite 
axial extension, + *>, Similarly, the plasnta density, p^ « Pp(r,c), 
current density, J s* J(r,a), and electric potential, « '^(r,a), are a 
superposition of primary (o) and secondary fields. 

The problem of secondary flows was first studied experimentally 
and theoretically for incompressible flow of liquids bet\ieen rotating 
cylinders (Taylor 19;i3) . A qualitative picture of the secondary flows 
in liquids between rotating cylinders is given in Fig, ?,2, which 
provides an impression of the complexity of these flows. In this case, 
radial and axial fluid motions have to be considered in addition to 


LL so that the primary fields (o) can be 




,uo 

a*s3fflttthal Qowponanta of fluid motiuu* ao tl\e Navle?»8tekea 

\ 

aquatiaua beeoma highly uonUuea?* Both espe^iroent CT^ayXoie 1923) aud 
hifuyoafeion theo^?y (tudovich 196?) tudieato that a oritioal Beyuolda 
numhat <\^, osiata Sqv \!M,q\\ tho psohlem Uaa a uuique aolutiou (Gouotte 
flow) for aroall Boyuolda uumhova I compared with < <\^), aud 
three«dim©UQlaual aolutioua of steady secondary flows after the loss 
of stability for I slightly larger than <X* S X)» As the. Reynolds 
number X ta f\srther increased » more and wore complicated types (n) of 
secondary flows occur as soon as X becomes larger than X^ (htfwoation) , 
where X are higher order eigen^valuss* Finally, for sufficiently 
large X, no laminar flow solutions exist and the rotating flow becomes 
ti\rbulent» The problem of secondary flows in ineotvpreasible flvtids has 
not been resolved analytically or numerically to date, 

Uve proposed plasm centrifxige model contains all major effects 
which t) enhance the rotation (electromagnetic forces, Uall effect) awl 
ii) reduce the rotation (viscous forces, bovindary layers, secondary 
flows) of the isotope mixture. Secondary flows occur because tbe 
oentrifvtge operates at high, awperaonlc speeds of rotation for which 
the Reynolds number X is probably larger than the critical value X^, 

The momentum and energy dissipated in the secondary flows represents a 
loss mechanism which reduces the speed of aslnxnthal plasma rotation and 
the efficiency of the centrifuge. For these reasons, the secondary 
flows have to be Included in the feasibility analysis of plasuta centri- 
fuges, As important is the determination of the critical Reynolds 
numbers X^ at which bifurcation occurs, which are important deaigu 
parameters for actual plaam centrifuges. 




...... -_-^= 



Analytical, llathod 

Itia two-Kiiina&a'lonal fields in the plasma eanti^ifuge with BQCon« 
dai*y flowa ace deacclbed by a nenlineac be«ndacy=valu© pceblem for the 
QQmpceaaible» maanetogaadynamle ©(luationo* 'Iliis problea\ ia Gon8ider= 
ably more Gomplex than that of the aeeondacy flowa of an incompreaalble 
liquid between rotating cylindero» ainoe the number of eoupled, 
nonlinear partial differential equations in compreaaible magnetogae- 
dynamics la much larger* ihe aolution of the nonlinear utagnetogaa^ 
dynamic equationa for the primary plasma field® ia first aougUt, and 
then the solution of the total plaan\a fields as a auperpoaition of 
primary and secondary fields* The secondary fields will be oKpanded 
in hiapunov’-ScUmidt aeries (Vaiidjorg and ‘iVenogUx» 1962) t e*g*, 

v(r,s) « S (X « X^)*^ V (r»a)» |X - X^| l» 
msl 

for the secondary velocity field* This statement for solution permits 
real aolutiona of the secondary flows for Ueynolds numbers X > 1^. 

Since |X - X^l « X in the vicinity of the onset of secondary flows, 

(X ~ Xy^)^ serves aa expansion parameter. Thus, in the vicinity of the 
onset of secondary flows, the notU.inear houndary"value proMem for the 
secondary fields can be treated by the method of successive approxima- 
tions. In this approach, the secondary fields in the m-th approx Im- 
tlon (m ^ I, 2, I, ... ) are described by coxipled, lixxear differential 
equationa, with eoefflcionto mxd source terms which depend only on the 
solutions of tixe lower approximations it I, 2, *,, m - I* In each 
approximation m ^ I to the coupled differential equations, the critical 
eigen-value X^ appears as an eigen-value, Tiuia, the higher 



U2 

nj»I>roisi\Batlovva alvo not only corVoctiona to tUo socondat'y flow flelclo 
but also to tUo eritical HoynoXds numboo Ttio aocondaty floldo 
havo to aatiafy boundary conditions at tho centrifuso valla i as well 
as periodicity conditions due to tho spatial poriodicity of tUo 
arrangemont of tho socondai'y flovs. Although this nxothod of solution 
is siffiplo in principle, tho actual integration of sovoval coupled, 
inhouKJSOttOQUQ difforontial equations is difficult, and the dogreo of 
cowplcnclty grows with each approximation m. 

It is proposed to analyao the plasma coutrifugo depicted in 
I'iguro A*l alovig the linos discussed above, in order to obtain 
solutions of tho plasma fields as a superposition of pvimry and 
secondary fio3.ds» l^ar tho mathematical details, it is referred to 
Section A, 3, which contains also physical extensions to the i^lasma 
centrifuge problem with secondary flows. 

7,3, i'Uoorotical l?ornvulation 

In tho theoretical description of the rotating plasma in the 
centrifuge, it is permitted tho inclusion of axisyawiotrlc secondary 
flows v(r,a) suporivwposcd on the main aniwvuthal plasma flow v^Q<r) so 
that tho velocity field of tho plasma is given by ^(r,») ® (v^,Cr,a), 
v^qCv) + VqCx'jb), Vj^(r,s)}. In view of tho extreme mtheawatical 
oomploxlty of tho av\alyois of secondary flows for comprosatblo centri- 
fuge flows, the assumptions of an isothoruvil partially-ionistcd plasm 
and small mgnotic Royinolds number arc owido first. .Accordingly, tho 
plasma tompocatv\ro is constant (T « 'f^) and tho induced magnetic fields 
can be neglected, » (0, 0, in this model. Tho proposed plasma 
centrifuge is described by tUo (isothermal) coaxprosslblo magnotogawd>n\a\wic 


equations and the conservation equation for the electric current 
density with Hall-effect (tox 0) , subject to the appropriate boundary 
conditions and periodicity conditions for the secondary flows. After 
normalization in accordance with the following substitutions. 


r/R r, z/R + z, 

^ (?.l) 

V’c'v 

where the characteristic reference values are defined by 


»c^Wc* 


(?.2) 


the dimensionless magnetogasdynamic equations for the mass density 
(p ), velocity (^), current density C?)» and electric potential (0) 
fields become 


» -Vpp-f- + -|vCv4)] + ^ ^ * (^*3) 

V*(p - 0 , C7./i) 

-wtJxS , (7.5) 

o o 


7 » ^ « 0 , ( J.6) 

where 

are Reynolds number and Hartmann number, respectively. Hie system 
(7.3)-(?,6) must satisfy the following boundary conditions. 






(7.7) 
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v„(r,a)rdr « 0, -» < a < + «* (.7.8) 

\ * 


** 2* (7.9) 

1)2 ^ 

where ^ end ^ dltnensioviless constants normalised by 

$ , respectively. 'Hie boundary conditions (.1.7) and (7.. 8) consider 
that the plasma does not slip at cylinder electrodes, and that there is 
no velocity flux tiiroush any transverse plane s *■ constant, respectively. 
Hio boundary conditions (7.9) specify that an electric potential 
difference '^2 “ is maintained across the perfectly conducting 
cylinder electrodes. Since the secondary flows are axlsymmetric and 
periodic along the cylinder axis (« •• axial wave number) the periodicity 
conditions for the secondary flows are 


Vj.^g(r, -E ±2it/a) - • 

<r, -a ±2v/ot) * “V„(r,n) . 

f* z 


(7.10) 


Hie formal perturbation theory is employed to analyse the 
beiiavior of the solutions in the vicinity of primary fields. All 
centrifuge ?(r,a) are therefore sought in the form 

« I^Cv) + f(r,a) , (7.11) 

Hie primary fields f (r) are first considered and it is supposed that 

o 

a small secondary field ?Cr,a) is superimposed on the initial steady- 
State field f^Ci^)* Substitution of the expression (7.11) into 









Equations (7. 3) “(7.9) yields a boundary-voluo problem for the initial 
steady state, and a linear eigen-value problem for the secondary flows 
by retaining only the first-order terms in the perturbation fields. 


Solutions for Primary Fields 


The following boundary-value problem for v^ (r) and 
describes the aero-order steady state of the plasma centrifuge; 


d fl d . „2 _ „2 ^'*‘0 

dr ^r dr ^ ^o0^^ " ^^08 dr ’ 


(7.. 12) 


r dr dr r dr oO^ ’ 


(f.l3) 


where 




(f.U) 




(-7.15) 


„2 _ „ 2 .., ,2 2 . 

H s H /(I + w T ) . 

1 


The remaining prinmry fields are consecutively determined by using the 


solutions for v _ and (|i : 
ou o 

991 

3„<r) - (1 + + v^e> • 


Jot ■ 


dr r ”0 “ * JoO • 


(7.16) 


(7.17) 


(7.18) 
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tl\e steady»statQ sttlutioRs o£ Equations (7»12)~(i?.,18) ars obtained In 
the forms 


v^gCr) « at 4* br*“^ 4- cr£nr , 


(J.19) 


.*2 #nx » • * .. 19.. _2 i 


- a(r‘/2) + bfur + c[rV2 - U )£itr - r /Al + d , <?,20) 

JL 




('. 21 ) 




or 


(T.22) 


f>.(r) » [2«4 tX r r ^esp(-/’'‘ r ^v^^diOdv + IJestpC/'’ r~^v\dr) , 

« li. 00 R, oo 


R. 




(^23) 


where 


a « -cUR^^hRj “rJ£iiR^)/<R^ -R^)1 » 


b - cCR^R^fiKRj/^^/C^ -rJ )1 


(7,24) 


c « -(vY^j)t(R^-Rj)/4 - R^R^en^(R2/Ri)/CR^-Rp + llf^MRj/E^) 


d « «*2-aR^/2 -bfnR, -ct(I^/2 - 2ir^)£i«R2 -R^/4] . 

JL 


Equations (?»20)«*(rf .23) exhibit clearly the effect of plasma cotat 




(?»19) on. the plasm fields. 


7.5. Sigeix-valuo Froblem for Secondary Fields 


Mnearisation of Equations (7,3)~(?.,9) according to Equation (7,ix) 
yields for the secovidary flows in the plasma the elgcn-vaiue problem: 


p (v *?u 4* v*?v ) 4- o(v^*Vv ) 
0 o o o o 


4- X“^E^^v 4- -| V(F*V)1 4> H‘X“^ j X 0 , 

.5 Z 
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!■h^! 
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1''1 

i \'..i 
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'in, 
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4l;n 

'■'v'1 


5 » •■ ■ , ■ 

. . ( ukWws**' 
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PqV’V + v*7p^ « 0 , 


(7.26) 


J -^ * 4 ' * 4 * 

•e -7^ + V >« e -WT j >< e , 


(7.27) 


V . J “ 0 


(V.28) 


»here 


v(r,z). 


< 2 < + ca 


(7.29) 


rRn 

/ * V (r,z)rdr «0, -oo<a< + «, 
Jo z 


(7 .30) 


4>(^^»z)-„o *“0» -»<2< + «. 

^ h,2 


(7 .31) 


B dcslgnaties the unit vector in the z~dlrection« Substitution of the 
2* 

periodicity conditions (7.10) for the secondary flows into Equations 
(,7,2S)-(7.28) produces additional periodicity conditions in z for the 
remaining perturbation fields; 


^r,6^^* -z ±2‘ir/a) = j^^g(r,z) , 


j (r, “Z ±2tr/a) « -j_(r,z) , 

Z A 


p(r, -z ±2ir/a)= p(r,z) , 


(J)(r, -z ±2v/o) =» <|>(r,z) . 


CV.32) 


In view of periodicity conditions (.7.10) and (7.32), the solutions 


for the linearized system are sought in the form 
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** \<^»«) • VjjCr) sin a« , 

j^^Q(r,n) “ 3^^0<s?) cos <*»» 3^(t,«> * 3^(r) sin ax , (7,33) 

p(r,«) « p(v) cos as* <^(r*s) “ |(r) cos a* * 

By substituting tho above trial solutions into tbo litioarlacd system 
(7,25)~(7,31) , tho problem is reduced to a sot of coupled ordinary 
difforontlal equations subject to homogeneous boundary conditions: 


(Dl)^ - “ ^'^“^1^0 ~ ^2^ ~ H^Cmt<-d| + Vg) -* Vj.3 


- i «»*\ + «»5,) • 

OM 

(DD^ - a^)VQ “ ■*■ ■*■ • 

JL 

(7.33) 

(Db* -* “ -Xap -f 1 + aVj^) , 

(7.36) 

(D^D “ (1 + ** + wTVp , 

(7.37) 

•** + av^ “ 0 . 

(7,38) 


vhere 



(7.39) 


(7.40) 


(7.41) 

♦<»>r-R2^2 ■ ® • 

(7.42) 

dJW^. * ® • 

(7.43) 



(7.U) 


9Xd 


XI 9 


d 

dr 




*j,(r) 5 2P„ »„8/r . »j<v) 5 vjj/r , tjW = , 


f^Cr) £ 


Tha perturbacton amplXtudas of the cxiri'enc density fields are 
determined by using the solutions for v^» Vq and 


3-Cr) « (1 + + Vfj •>* WTV ) , (?’A5) 

VJ i 

JgCr) ** (1 + * <“.A6) 

3 C^*) ct| . (?«A7) 


Equations with the corresponding boundary 

eonditions (?.39)-(7.A3) constitute an eigen-value problem for the 
Eeynolds number ,\, For given system constants 2 » h» 
operating parameters (H, 'tg), the eigen-valv^e problent (?.3d)-(7 .A3) 

has a sequence of positive eigexi-vaXues fnt' any axial vtave 

number «» Among these eigen-values the snvallest one, ” 

Q^<» ^ certain value of is the critical Reynolds number, 

at which the secondary flows first set in. 

The similar eigenvalv\e problems as sbo\<n here have frequently 
occurred in the analysts of axisytimetric neutral stability for either a 
hydrodynamic Gouec^e flow or an IffllD Gouette flow. Indeed even In these 
simpler cases, the problems have never been completely solved by 
analytical methods. Either case has been analysed using soma 
techniques (Ghandrasehhar 1961, Ovchinnikova and ludovich 1966) of an 
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expansion in ozthogonal functions, a variational method and Green's 
function solution or direct numerical computations by employing the 
Galerkin method (Kurzweg 1963), the Runge-Kutta method (Harris and 
Reid 1964), successive approxiniatlons (Sparrow et al, 1964) or a 
finite-difference technique (Chang and Sartory 1965) . 

The solutions of the eigenvalue problem (?,34)-(7,43) correspond 
to the first approximation (m=l) of the Liapunov-Schmidt series 
expansion (Section 7.2). In a similar way, the eigenvalue problems 
for the higher order expansion fields v (m>2) can be formulated to 
obtain improved series solutions and eigenvalues X for the secondary 
fields. Along these lines, a quantitative theory of. the plasma 
centrifuge with secondary flows could be developed with an accuracy 
corresponding to the second approximation (m»2) of the Liapunov-Schmidt 
expansion. 

In the above theoretical formulation of the plasma centrifuge 
problem, certain physical effects have been neglected in order to 
reduce the number of equations and to reduce the formalism. In the 
actual research, induced magnetic fields should be taken Into consid- 
eration so that the results are applicable for arbitrary magnetic 
Reynolds numbers, R ® large thermal 

conductivity of plasmas, the heavy particle temperature is quasi- 

homogeneous over most of the centrifuge space R^ < r < i«c.» 

temperature drops exist in the vicinity of the electrode walls. Thermal 
energy transport and dissipation should be taken into account If this is 
mathematically feasible. The theoretical approach proposed would make it 
also possible to solve the dynamics of secondary flows in ordinary, 
nonconducting gases, another still unsolved problem. 
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aiAPTER VIII 

COLUSIONLESS PLASM CENTRIFUGE 

In this section, the problem of the collislonless plasma 
centrifuge is presented within the framework of the (steady-state) 
Vlasov-tfesrs'^ell equations. A one-dimensional configuration of a 
multicomponent rarefied plasma in an infinitely long cylinder aligned 
parallel to an external axial magnetic field is considered. In the 
charge-neutral approximation! a solution to the self-consistent 
VlasovHaxijell equations is obtained by assuming Masr^ellian distribu- 
tions for the plasma particles. 

8.1. Mdel 

As a model for a collisionless plasma centrifuge an electrically 
insulating cylinder R^ containing a mixture of isotope ions and 
electrons of toasses m^, i » l!2,..,e is chosen* The plasma state is 
produced by s\dltching on a strong axial magnetic field B(t) ->• - 1 

Tesla so that the associated induced electric field E(t) , which is in 
the asimuthal direction! breaks down the isotope mixture simultaneously. 
The ions and electrons are accelerated by the Induced electric field so 
that they obtain mean mass velocities in the azimuthal direction. 
The resulting centrifugal forces (in a system of reference moving with 
the particles) distribute the Isotope ions radially in accordance with 
their different masses m^* The subject of the consideration is the 
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steady-3 tate composition of the isotope itiixture after the laagnetlc 
field B(t) has reached its "plateau” value (Fig. 8,1). 



Fig. 8.1. Model of collisionless plasma centrifuge. 

This centrifuge has to be operated at low particle densities 
(collisionless system without significant viscous losses at r = R^) , 
but at still high enough pressures to avoid a compression of the plasma 
Immediately after breakdown by the Lorentz force JqB^ which is 
directed radially inwards (theta pinch effect). A pinch contraction of 
the plasma with a radial velocity ^ would produce even higher azimuthal 
particle velocities through the azimuthal V x B field but would render 
a controlled extraction of the isotope ions difficult within the typical 
pinch times At - 10 sec. Another reason for avoiding 0-pinch 


CQwdifeions is the sUock Uesting of tUa iona to tempovaturea 10^ ®K. 
At such high ion ten^pcvaturoa » the tendom thei^mal forces would dominate 
the directed centrifugal forces and make an efficient Isotope separation 
impossible , 

Accordingly, in an e;?perlTOent the centrifuge would have to be 
operated under conditions where the radial Lorenta force ia small 
compared to the pressure gradient and the electron temperature 
generated in the induced electrical breakdown is large compared to the 
ion temperature, 10 “K, Nonisotherwal plasmas are 

readily realised at low filling pressure T^) since only a 

fraction - m M, of energy is lost by an electron e in collision with an 
ion i« Since ideal collision and lossless plasms do not exist, only a 
quasi-equilibrium can be reached in the isotope ion mixture after B(t) 
has risen to which shovild last for about At ' 10 sec, and permit 
extraction of the separated isotope species due to the long thermal 
relaxation tia\e of the electrons. As a technical application, this 
collisionless plasma centrifuge would have to be operated under 
continuously repeated Induction pulses h(t) '>* in order to separate 
a significant amount of Isotopes in a reasonable time. 

8,2. Boundary-value Problem for Vlasov-Maswell Equations 

From the theoretical point of view, one is interested in analyssing 

the quasi-equilibrium state in the final external magnetic field 

and the self-consistent electric, E » and magnetic, B » F x A 

fields. For a quasi-inf initely long centrifuge, all fields are functions 

of V only, i.e., E « -e„3({!/3r and B« o a(rA)/rar. In the one-dimensional 

V * 
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steady state, the distribution functions f^(^,r) of the particles of 
type i * l,2,.,,e are described by the coupled Vlasov-Haxwell 
equations (Section 2.2) : 


1 a OA vf 3f. 

y — i 4- [_i Y i „ _JL .M. + — i. 

8 r ^ ''e r 3r^*^ ^ in. 9r r ^ 3V 

i i r 


.. r— i. Y i _L.fj-A) 4- - -■ -] — i => 0 . 
r r 3r'- '' r ^ 3 Vq ’ 


( 8 . 1 ) 


where 


X ® 


(S.2) 


13 3^.,^ _-l 


Ofl 


■— ^(r e £ e * /// f * dV , 
r 3r 3r 0 , i i ’ 

1 - “ 


(9.3) 


relate the vector potential A(r) and the scalar potential iji (r) to its 
current density and space charge sources, respectively (V , V„, V 
are the cylindrical coordinates of the particles 1 in the velocity 
space) . 

The equations of motion for each particle species 1, 




dV 
r 

dt 


-e. 


f ^ "l’^67 


3r 


(rA) + “ 


1 a • 

m. -TT" “ ~e ,V — -^(rA) 
i dt i r r 3r 




V V. 
r 6 




dt 


(8.4) 


have the three integrals of motion: 




125 


Hi - + Vq + V^) + e^<|. , 

(8.5) 

Pi 5 - m^Vgr + ej^Ar , 

(8.6) 


(8.7) 

The Hamiltonian the generalized angular momentum 

H P^g and the 


generalized axial momentum are constants of motion for each 
isotope i. Knowing the complete set of integrals of motion, the 
general solutions of equations (8,1) are arbitrary functionals of 
and does not enter explicitly due to the absence of the mean 

axial velocity) : 


expC-a^Hj^ - . 


( 8 . 8 ) 


The arbitrary constants and 3^ are to be determined by 

Max^^ellian boundary conditions at r “ 


f^C^,r) 


„ t ^ 3/2 

Vo.trU'P ' ^ 


r=R "oi'2TrkT. 

o i 


exp{-m^[vj + CV^oi^^ > 


(8.9) 


S T^, i « l,2,,..i»e; T^. s 1 


e 


Note that n^^ and are the particle densities and mean (azimuthal) 
mass velocities (rotation) of the i-particles at r = R^. Since 

3 

T » T due to electron heating and T - 10 “K is of the order of 
e 0 ”0 

the wall temperature, a distinction between temperatures 1 1 and J_ to 
is unnecessary. 
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Substitution of Equation (8.9) into Equation (8,8) yields 


« 1/kT^ , 


• 


( 8 . 10 ) 


‘Bi. 3/2 

^Oi '* \l^2irE7^ exp[-(— y 

1 


- Vo + 


M. = V ./R , = (}i(R ), A s A(R ) 

i oi o ^0 o'* o o 


(8.11) 


Xhe particle density and aaimuthal mean velocity of each isotope are 


given by 


OQ 

n.(r) « /// f.d^ 

* ... w ^ 


2Ttk’T 


2 2 
in^ta^r 


*ol< m > 0=q.t(-^ - + Cj^<«jrA)/kT^] <8.12) 


’iW - ^ ■ V • 


Substitution of Equations (8.12) and (8,13) into Equations (8,2) and 


(8.3) yields a boundary-value problem for coupled nonlinear ordinary 


differential equations for the self-consistent potentials: 


^t|^(rA)]--M„r 2 ■ 


i^rf) I Oj^n^(r) , 
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* VA > 


ill 


dl*’ e 2tte R f * 

0 0 i 
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fl "O. 


and 


r«o 

»* S 2w J tt. i*dv 

i '^0 i 


(8.16) 

(8.17) 

(8.18) 
(8.19) 


(8.a0) 


is the total nurabet o£ 1 particles per milt length. The boundary 
conditions f8.16)~(8,19) can be derived directly from Bquations (8.14) 
and B^uatlon (.8,15) by integration > respectively. The boundary 
conditions (.8.16) and (.8.17) consider that the requirements of 
syiMsietry make the self-consistent magnetic field parallel to the 
cylinder astis at points inside the centrifuge and that it approaches 
aero at points outside the centrifuge, the boundary conditions (8.18) 
and (8,19) take into consideration that the radial electric field is 
due to the space charges. 


8.3. ()uasi-neui.i«l sslutions 


In the limiting case of charge neutrality (X e^n^j^ » 0), there are 

i 

no internal space charges so that the electric field vanishes and only 
a magnetic field exists. In this case, by Equation (8.12), 


2tkt 

g A ^ 

t ^ ®^i 


1,3/2 

A) expK — 




(8,21) 
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from which 




a -( const, s b « const , 

StvIcT, 

1 •*“ «* “i 


( 8 . 22 ) 


Than, Equation (8.14) becomes 


^tf^CrA)] . -U^r exp(b/ + arA) S “l ' 


2TrkT, 


With the change of variables 


s br^ + arA , 




(8.24) 


where 


9 , 2ffkT. 3/2 

- 1 >-0^ I “i 

f! . i i ■ 


(8.25) 


Equation (8 ,23) can be transformed to 


dt 


(8.26) 


ihe general solution of Equation (8.26) has the form 

4.,_sare^‘‘_ 


(8.27) 


tdiere G, \ and y determined by Equation (8,25) and the appropriate 
boundary conditions: 






Y ** - Q)/(p + Q) 


(8,28) 


X « 




t) , 


vliero 


I Wi * 

Q 2 iUB^ + + 2 Uq S 


(8.29) 


(8.30) 


Is the magnotte field vhicli is necessary for yavticlo i to rotate 
with tlio eyclotrou froquoucy It is noted tluat 


(2ir)”^ £ *^i**'i^i “ ^ total current per unit length, and that 

£ tt j.kT. e £ n whore V* is a mean square velocity of each i. By 

l^oiij^oiir r 

the above expressions, tho quasi-noutrol solutions are finally given by 


(}sCr) •• 0, V^(r) “ Mj^r, J^(r) « 


* ”oi '3 




Jta 


Xr 


.2 




\r2 

B(r) - -B* + Q . 

1+YO 


(8.31) 


It turns out, however, that tho radial distribution of the 
isotope ions depends sensitively on tho electric potential 't'(r) since 
Oj't'/UT > 1 for the ions i j* o(kT ~ 75*^10"^ volt for T » 10^ *K), 

JL vJ VJ %J 

2 

i.e,, it Is not determined solely by tho centrifugal forces ^’oi 
Por this reason, tho theory should bo extended to include space charge 
effects and to obtain tho exact potential distribution ^JiCr) . In view 
of the complexity of tho underlying i\onlinear (transcendental) dlt I'ev- 
ontlal equations, this extension has probably to be eavrted tlrough 
vvflthln tbe frame of an appropriate pertnvbat Ion theorv. 


APPENDIX: PREVIOUS RESEARCH ON ISOTOPE SEPARATION 

a 

For isotope separation* various methods have been suggested such 
as chemical methods (Urey 1939), diffusion methods (Furry et al. 1939), 
electromagnetic methods of mass spectrometry (Smith et al. 1947), and 
mechanical centrifuge methods (Humphreys 1939 j Cohen 1951). In electro- 
magnetic methods, ions are moving in different orbits and are deflected 
in a magnetic field according to their different charge to mass ratio. 
However, because of difficulties in producing Intense ion beams and 
neutralization of ions by electrons, this method was not widely used in 
industry. In fact, it is used only for laboratory purpose for producing 
limited amounts of pure Isotopes owing to its lilgh resolution. For 
separation of isotopes with low mass, chemical methods are more effec- 
tive than electromagnetic methods, Nlien large quantities of pure 
materials are required, it is Industrial practice to use gaseous-diffusion 
separation systems. 

In recent years, other effective methods for isotope separation have 
been studied. In particular, plasma methods are promising for high pre- 
cision technology and in new technological developments. Plasma separa- 
tors can be designed to operate on the basis of plasma streams (Becker 
nozzle) and plasma rotation (plasma centrifuge). The Becker nozzle is 
being developed as a commercial isotope separation device. In principle, 
this device expands the isotope mixture through a supersonic nozzle and 
along a curved wall so that extremely large centrifugal forces result 
which separate the heavy isotope from the light (Becker et al. 1955). 

An extension of this principle is applied in the jet scheme in which two 
or more opposing supersonic nozzle flows deflect each other so that 
centrifugal forces occur again as a result of stream line curvature 
(Campai'gue 1970} Becker et al. 1973). 
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Although the plasma centrifuge concept was first proposed by Slepian 
(1956), a major research effort does not appear to exist in this field in 
the United States. There is a significant research program on plasma 
centrifuges in the U.S.S.R. (Berezov et al. 1976; Belorusov et al. 1976), 
which is classified. The interest of smaller countries, such as Sweden 
(Bonnevier 1966, 1971; Lehnert 1970, 1973), Japan (Okada et al. 1973), 
and Australia (George and Kane 1972; James and Simpson 1974, 1976), in 
plasma centrifuges appears to be due to the low cost of this type of 
separation device. From the theoretical point of view, the basic mecha- 
nism for plasma rotation by means of crossed electric and magnetic fields 
and Lorentz forces in rarefied and dense plasmas is vinderstood qualita- 
tively (Anderson et al. 1958; Gordeev 1959, 1961; Kessey 1964; Hanson 
and Cohen 1970; Vrba 1971; Witalis 1974; Ban and Seklguchi 1976; Marlier 
1977). 

Proposed plasma centrifuges employ either low-density collisionless 
plasmas or high-density collision-dominated plasmas as working fluids. 
Experimental evidence on isotope separation in plasma centrifuges has 
been reported for both cases (Bonnevier 1971; James and Simpson 1974, 
1976; Heller and Simon 1974; Berezov et al, 1976), Exact solutions for 
collisionless centrifuge plasmas are not known, which require evaluation 
of the self-consistent electromagnetic field interactions (Komarov and 
Fadeev 1962; Watson 1956). The disadvantages of collislonless centri- 
fuges are relatively large electric power dissipation to produce high 
degrees of ionization of the isotope mixture and the small amounts of 
isotopes they permit to separate. On the other hand, collislonless 
centrifuges have minimum velocity losses at the walls due to the absence 
of ordinary hydrodynamic boundary layers. In high-density centrifuges. 



only a small fraction of isotope such as cesium with a low ioniaation 
energy) has to be ionised to produce a partially ionised plasma state. 
The horents force due to the interaction of the cvurrent density and 
magnetic field sets not only the charged plasma components but also the 
neutral plasma components in rotation* which are coupled through the 
intorcomponent friction forces* The collision-dominated centrifuge 
excels through relatively low energy dissipation and largo isotope 
densities. The velocity losses occurring in the viscous boundary layers 
at the walls are* however* of some disadvantage. 
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